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FOREWORD 

This report presents results of work performed 
by Lockheed's Huntsville Research & Engineering 
Center for the NASA-Mar shall Space Flight Center 
under Contract NAS8-28170, "Methods for Structural 
Design at Elevated Temperature". 

The NASA-MSFC technical monitor for this 
contract is Mr. John E. Key of the Analytical Me- 
chanics Division, Astronautics Laboratory. 
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SUMMARY 

The factor of safety used to compute the ultimate load from the limit 
load is a factor to account for uncertainties both in the strength and stiffness 
of the structure and in the applied loads. This uncertainty is twofold; first, 
in the physical reality and second, in the mathematical description. 

The objective of this study is to define a procedure which can be used to 
design elevated temperature structures. The desired goal is to have the same 
confidence in the structural integrity at elevated temperature as the factor of 
safety gives on mechanical loads at room temperature. 

Methods of design and analysis for creep, creep rupture, and creep 
buckling are presented in Section 2. Example problems are included to 
illustrate the analytical methods. Creep data for some common structural 
materials are presented in Appendix A. Appendix B is description, user's 
manual, and listing for the Creep Analysis Program developed on the contract. 
The program predicts time to a given creep or to creep rupture for a material 
subjected to a specified stress-temperature-time spectrum. 

Fatigue at elevated temperature is discussed in Section 3. Methods of 
analysis for high stress-low cycle fatigue, fatigue below the creep range, and 
fatigue in the creep range are discussed. The interaction of thermal fatigue 
and mechanical loads is considered, and a detailed approach to fatigue analysis 
is given for structures operating below the creep range. 

Section 4 is a brief discussion of structural analysis at elevated tem- 
perature. Limitations of linear, elastic analyses and the desirability of 
developing efficient nonlinear analytical tools are pointed out. 
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Design methods for elevated temperature structures are recommended 
in Section 5. Both ultimate load failure (fracture) and functional failure from 
excessive permanent deformations are considered as design criteria. The 
recommended approach to elevated temperature design consists of applying 
a factor of safety to the mechanical loads, a life factor to the service life, 
and a factor to the design heating rates or temperatures. The rationale for 
applying the factors is explained based on reliability principles. The design 
load, temperature, and life factors assure that the design conditions (load, 
temperature and time) have a safe margin when compared to the allowable 
(failure or deformation) envelope. 
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Section 1 
INTRODUCTION 

In their simplest form, the criteria for designing room temperature 
structures are: 

1. The limit load is defined as the maximum load expected to act 
on a structure; 

2. The limit load is multiplied by a factor of safety to obtain the 
ultimate load ; 

3. The structure must possess adequate strength and stiffness to 
withstand limit load without yielding; and 

4. The structure must withstand ultimate load without failure. 

The factor of safety used to compute the ultimate load from the limit 
load is a factor to account for uncertainties both in the strength and stiff- 
ness of the structure and in the applied loads. This uncertainty is twofold; 
first, in the physical reality and second, in the mathematical description. 
Strength will vary with material properties, manufacturing tolerances, and 
fabrication techniques. The applied loads vary with trajectory perturbations, 
wind gusts, maneuver loads, atmospheric conditions, and numerous other 
random phenomena encountered during service. 

Due to scatter in measured material properties and the randomness of 
the loads, which are highly probabilistic in magnitude and frequency of occur- 
rence, it has been logically hypothesized that the factors of safety should be 
determined from a statistical analysis of the basic factors which affect the 
strength and the loads. Unfortunately, a statistical analysis is usually not 
possible in an early design stage and also it is too time-consuming to be 
performed in every step of a structural design. Therefore, design criteria 
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are desirable which convert the probabilistic properties of the design variables 
into deterministic ones. This leads to the factor of safety, which has been used 
as an effective design tool for a long time. 

Although the factor of safety approach to designing reliable structures 
has served well, the high temperatures to which space shuttle structures will 
be exposed require a generalization of this approach to a wider range of appli- 
cation. At elevated temperatures the uncertainties of the design variables are 
greater and call for larger factors of safety as compared to room temperature. 
In addition, the factors of safety should be functions of the design variables 
rather than a constant. The latter aspect is emphasized by the large percentage 
of vehicle weight incorporated in elevated temperature structures. Excessive 
conservatism as a result of an insensitive factor of safety would buy structural 
reliability at the expense of high weight penalties. 

The objective of this study is to define a procedure which can be used to 
design elevated temperature structures. The desired goal is to have the same 
confidence in the structural integrity at elevated temperature as the factor of 
safety gives on mechanical loads at room temperature. 

The problem of design of elevated temperature structures is to attain a 
safe, reliable structure for combined load-temperature-time environments. 
Elevated temperature design differs from "room temperature" design only when 
the combined environments and material properties interact to produce changes 
in the strength or deformation characteristics of the structure as the load or 
temperature cycles are repeated. That is, the time becomes a major param- 
eter in the design of elevated temperature structures. 

The following sections of this report discuss methods for structural de- 
sign and analysis which are peculiar to elevated temperature structures: creep, 
creep rupture, thermal and mechanical fatigue, and thermal stresses. The 
interaction of the thermal and mechanical loads to produce accumulated damage 
and permanent deformations with service lifetime are considered. An approach 
to design of elevated temperatures is proposed. 
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Section 2 
CREEP 

This section presents a summary of the phenomenon of creep and its 
effects on the design of a structural component at elevated temperatures. 

Creep is one of the major drivers in elevated temperature design, especially 
in flight structures where permanent deformations are undesirable from an 
aerodynamic standpoint. This section covers the mathematical theories of 
creep from the aspect of material laws and description of the creep behavior 
of structural components. Questions of creep rupture and creep buckling are 
discussed. Finally, a computer program is described which implements three 
popular creep theories used in studying the effects of a time -temperature -load 
history on a structural component. 

2.1 CREEP LAWS 


The mechanisms of deformation and failure of metallic materials at 
elevated temperature, that is, at temperatures above approximately one -third 
of the melting temperature measured on an absolute scale, are traditionally 
studied on two levels, the crystal lattice and the continuum (Ref. 2-1). 


On the crystal lattice level there are two softening mechanisms, cross 
slip and the temperature -sensitive dislocation climb. If the stress is kept 
constant the hardening effect decreases while the softening increases rapidly 
as the temperature is raised. The strain rate, de/dt, can be expressed in 
terms of the stress, a, as 


de 

dt 


~ n 
C a e 


_Q 

RT 


( 2 . 1 ) 


where Q is the activation energy for self diffusion, R is the gas constant, and 
T is the absolute temperature. There are several empirical methods to relate 
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a, T and time, t, at rupture or at a given strain. A frequently used parameter 
is the one by Larson and Miller (Ref. 2-2), where 

P(o) = T (C + log t) . (2.2) 

A list of parameters is given on page 200 of Ref. 2-3, and methods are described 
for extrapolating creep and creep rupture data. Typically the data are in the 
form of "master curves" (Fig. 2-1) from which individual creep curves of creep 
strain versus time at constant temperature with parametric variation of stress 
(Fig.2-2a) or creep strain versus time at constant stress with parametric varia- 
tion of temperature (Fig.2-2b) can be constructed. In Appendix A the sources 
of material data and creep curves for a number of metallic materials are pre- 
sented. 

The material behavior shown on Fig. 2-2 can be described approximately 
by 

e • (tP *f (%)' * ■ 

o 

as shown on Fig. 2-3. 


e 



Fig. 2-3 - Approximation of Creep Curve 
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The term 


■ (t r 


(2.4) 


(Norton) describes the plastic strain with strain hardening, while 


de 


dt 


(tf 


(2.5) 


represents the viscous flow under constant stress. The elastic strain 


a 


( 2 . 6 ) 


has been neglected in Eq. (2.3). 

Before the mathematical implications of Eq. (2.3) are pursued any further 
it may be of interest to explore alternate formulations. Generally all strain 
theories can be written as 

e =f 1 (a)f 2 (t)f 3 (T) (2.7) 

in a multiplicative form of functions of stress, <r, time, t, and temperature, T. 
Particular "rules" are the time hardening, 

d€ d f 2<*> 

af ■ f i<»> -ar- £ 3 (T > < 2 - 8 > 

strain hardening, 

ff- = g x (a) g 2 (e) g 3 (T) (2.9) 

or a combined theory as 
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= C a a e P t r (2.10) 

Other material laws are expressed as 

CT P 

C — (Nadai) (2.11) 


(Cj q e‘ qt + C 2 ) (J a (Marin). (2.12) 

The material law of Eq. (2.3) has advantages for steady state problems 
(secondary creep) while for transient problems (primary creep) the material 
laws of Eqs. (2.7) through (2.12) have to be considered. 


and, 


de 

dt 


de 

dt 


To complete the catalog of material laws Eq. (2.3) is generalized to loading 
and unloading problems, including recovery effects (Ref. 2-4). During loading 


de 

dt 


_d_ 

dt 


l* ♦(*)“>(*)■ 


(2.13) 


while during unloading 


de 

dt 


_d_ 

dt 


l* * m + m n 


(2.14) 


This creep law is illustrated on Fig. 2-4. 


2.2 CREEP RUPTURE LAWS 


Creep rupture is a failure condition which is associated with a particular 
stress temperature and rupture time. Since Hencky 1 s notable paper of 1925 
(Ref. 2-5) a number of theories to predict rupture times have appeared (Refs. 
2-6 through 2-9). These papers are briefly summarized on Fig. 2-5, on which 
the following notation is used: 
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Name 

Material Law 

Time to 

Ductile Rupture 

Time to 
Brittle Rupture 

1. 

Hencky 

1925 

<Oy = dy 0 + 3KE 

t R « — In 

yo ^t> ~ G'ya 


2. 

Odqvist 

1933 

G* y a ^ y o t 3 ice 
+ Ke 

3 fc - <^o ~ K 

to = '•» 

~ <^yo 


3. 

Robinson 

1952 

'f (T e ,-tO 

Larson-Miller 
Ma n s on - Haf e rd 

t R 

5 f - i 

o T * 

4. 

Hoff 

1953 


v i (ft)' 


Hoff with 

primary 

creep 

’ f <o \ / \ 

€*/“ + 7 

y y 



6. 

Kachanov 

1958 

C^V- — 

R * r+l 

t K 

S(£' 

O 

j dt - 


Fig. 2-5 - Creep Rupture Laws 
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yield stress 

initial yield stress 

"viscosity" 

initial stress (at t=0) 

material constants for incipient creep 

material constants for steady state creep 

material constants for creep rupture 

rupture time. 


Robinson's theory is very useful in predicting rupture times for complicated 
temperature-stress histories since it applies the law of linear cumulative 
damage to creep rupture. To compute the expended lifetime to rupture the 
material data from a "master curve" such as the Larson-Miller curve are 
evaluated and the contributions of the individual phases of the temperature- 
stress history are summed. 


Kachanov's theory can be combined with the theories of ductile rupture 
in lines (4) and (5) of Fig. 2-5 so that a prediction of both ductile and brittle 
rupture is possible. The terms ductile and brittle refer to the amount of cross- 
sectional necking. When the cross-section shrinks to zero at rupture, it is 
called ductile. Kachanov's theory can be used to calculate rupture times for 
multiaxial states of stress in which either the maximum tensile stress or, as 
may be the case in some materials, the equivalent stress a is used. In the 
case of non -homogeneous stresses the travel of a rupture surface through a 
structural member can be treated using Kachanov's theory. 


Another form of presenting Kachanov's theory is as follows 


where 


dD 

dt 


C 



n 


D 



(2.15) 


( 2 . 16 ) 
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is the "damage factor" which approaches 1 when the reduced cross-sectional 
area approaches zero. 

2.3 EVALUATION OF MATERIAL CONSTANTS 

The material law given by Eq. (2.3) is written in a slightly different form 


e = Cj tf m + C 2 cr n t (2.17) 

The material constants Cj and m are associated with the incipient (or intercept) 
strain e . and C-, and n with the steady state creep strain e (Fig. 2-6) 

1 u C 



Fig. 2-6 - Simplified Creep Law 


The intercept strain data are presented in the form of Fig. 2-7, where 
the actual curve is approximately a straight line with the slope m. Then 
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log a 



Fig. 2-8 - Master Creep Curves 


and at stress level 



To find n and the relation 



is used so that 


( 2 . 22 ) 


(2.23) 


n = In (e j/e 2 )/fn (ctj/o^) 


(2.24) 


and 



(2.25) 
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2.4 LIFETIME CALCULATION IN A RANDOM TEMPERATURE FIELD 


In the following the creep law is considered in the form (Ref. 2-10) 


Sr ■ K °“ 


(2.26) 


where for moderate temperature variations 


K = c e 


P T(t) 


(2.27) 


This is the "viscosity factor" with c and (3 being material constants. 


When the temperature is subject to temperature fluctuations 


T(t) =T m (t) + 0(t) 


(2.28) 


where is a mean temperature and 9 is the random fluctuation, then the 

viscosity constant becomes sensitive to these fluctuations 


K = c e p ' x m 


P(T m (t) + e (t)) = pe(t) 

m 


(2.29) 


where 


K = c e P T m^ (2.30) 

m 

If the probability density function p(0) of the temperature fluctuation is known, 
the expected (mean) value of the viscosity constant is 

+00 

K = K m / e (39(t) p(0)de = K m Mt) (2.31) 

-00 
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For a normal distribution 


M(t) = e 2 (2.32) 

2 

where T is the variance of the distribution. 

This formulation can be used to find the reduction of lifetime of structures 
subjected to temperature fluctuations in the high temperature range. 

Consider a bar under constant load P. The necking of the bar expressed 
as 

t ^ (2.33) 

o 

where A q is the original cross-sectional area, is a measure of the lifetime 

1 >;_>° • (2.34) 

It is assumed that the life is up when t, = 0. 

The natural strain is 


and the strain rate 


while the stress is 


where ct q 
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Using Norton's law in the form of Eq. (2.26) the following differential equation 
is arrived at 


t + c a. 


e 0T(t) ; l-n 


which leads to the integrals 


C (t) 

J ; n_1 d t, 

U°) 


t 


-c 



e 0T(t) 


o 


dt 


(2.38) 


(2.39) 


Taking expectations on both sides of Eq. (2.36) and setting the expected value of 
t, at lifetime to be zero the following result is obtained 


t 


L 



M(t) K m (t) dt = 


o 


1 

n 

nCT o 


(2.40) 


For stationary conditions the lifetime can be explicitly expressed from Eq. (2.40) 
by 


*1 = 


1 


n 

n a u K 

o r m 


(2.41) 


2.5 VARIATIONS IN STRESS AND TEMPERATURE 


It is instructive to qualitatively examine the behavior of a constrained 
bar with various material laws involving creep effects, when it is subjected 
to temperature cycles (Refs. 2-11 and 2-12). Figure 2-9 is an extension of the 
case of elastic -plastic material with constant yield stress and of elastic -linear 
work hardening with various types of yield strength (temperature dependent, 
unstable, etc.) involving the accumulation of stress (residual stress). The 
establishment of various cyclic patterns according to heating and cooling rates 
is evident, Fig. 2-10, 
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(a) Viscoelastic-Viscoplastic 


6 = 



a 



a - a (T) 
+ - 


(b) Slow Heating and Slow Cooling 


. a a 

6 = E + cY 


Fig. 2-9 - Thermo-Cycling of a Clamped Bar 
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Fig. 2-10 - Viscoelastic-Plasti 
Cooling Rates 
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The case of variable stress at constant temperatures is illustrated on 
Fig. 2-11. The accumulation of creep strain when calculated with some of 
the theories discussed earlier, Eqs. (2.7 through 2.12) is depicted. In Ref. 
2-3 a number of stress histories and their effect on the strain accumulation 
are discussed. 

2.6 MATHEMATICAL THEORIES FOR ANALYSIS OF STRUCTURAL COM- 
PONENTS SUBJECTED TO CREEP 


Consider the tensor of deviator stresses S„, where i, j = 1,2,3 refer 
to a rectangular Cartesian coordinate system. For isotropic materials an 
equivalent stress is defined by 


a 2 = jS.. S.. 
e 2 ij ij 


(2.42) 


The summation convention of repeated indices is used. Similarly an equivalent 

strain e is 
e 


2 2 

e = 

e 3 ij i] 


and an equivalent strain rate v g is 


where 


2 


v 

e 



v. . 
ij 


v. . 
U 



(2.43) 


(2.44) 


(2.45) 


It is assumed that the material law of Eq. (2.3) is used. Then 



(2.46) 


and 
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The corresponding tensor relations are 


and 



(2.47) 


(2.48) 


(2.49) 


The bending of a beam following the material law of Eq. (2.47) can be described 
in a similar fashion as an elastic beam. This has been called the elastic 
analog (Ref. 2-4). The strain rate is 


v 


9 3 w 

~ z 2 

8 x 9 t 


* n 

-z w 


Rewriting the material law as 


v 



sgn a 


(2.50) 


(2.51) 


and solving for a gives 


a = a c |v| n sgn v 


(2.52) 


Let a bending moment M and an axial force N act on the cross section (Fig. 
2-12). Then the stress resultants in terms of the stress are 


e 2 +d 


/ 


ffdA = N 


(2.53) 


-©j + d 


2-19 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D306579 




Fig. 2-12 - Beam Cross Section 


e 2 + d 



a (z-d) dA 


-e j + d 


M 


which can be rewritten with Eqs. (2.50) and (2.51) as 


a 

c 


v 


1 _ 

n e 

sgn v S n 


N 



1 _ 

11 T 

sgn v I n 


M + dN 
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The area integrals are 


*2 + d . 1 


-e j + d 


n 


sgn z b(z) dz 


e 2 + d 


I = 
n 


1+i 

n 


b(z) dz 


-e ^ + d 


(2.57) 


(2.58) 


On Fig. 2-13 the area integrals for a few cross sections are given. 

The material law of Eqs. (2.46) and (2.47) also lends itself to deriving 
a beam theory, as well as plate and shell theories. A beam equation similar 
to Bernoulli's equation for elastic beams can be visualized from Eq. (2.56). 

For example a cantilever beam with a tip load P has a bending moment 

M(x) = P x (2.59) 


if the x-axis runs from the tip. 


From Eq. (2.56) it follows that 


v = 


. ii 

-w 


O' I 
' c n 


. n 


n 

x 


(2.60) 


It is a simple matter to derive beam solutions for a variety of boundary condi- 
tions (Fig. 2-14). 


Using the foregoing the lifetime of a cantilever beam with a tip load, 
subjected to the temperature fluctuation of Eq. (2.28) would be computed from 
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Fig. 2-13 - Area Integrals for Beam Cross Sections for 
Elastic Analog 
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Fig. 2-14 - Deflection Rate for Beams with 
Various Boundary Conditions 
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^ 


\ n T 2 

- n K t 

. I n+Z ^ m 
n ' 


+ (2i_\ r 

3EI + l I„ / 

when a deflection limitation is given. 

Another application is the development of design charts for combined 
loads. The largest effective stress for a non-linear structure subjected to 
r loads L^, L^, . . . , acting simultaneously, 


( 2 . 61 ) 


^ * f n < L 1- L 2- 


• L r> 


(2.62) 


or with 


L. 

a- = y~ i = 2, 3, 
x 


. r 


(2.63) 


o ^ Lj f^ ( 1 , 0^2* otj , • • • » oc r ) 


(2.64) 


can be compared to the effective stress for an equivalent linear structure 


f j (l»a 2 » a 3» * * * ,a 4 ) 


(2.65) 


by introducing (Ref. 2-13) 


0 



( 2 . 66 ) 


If tables or charts of 0 are available, then by simply computing cr with standard 
handbooks method and looking up 0 the equivalent stress of the nonlinear struc- 
ture is 


<>>= 0 


(2.67) 
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For a beam subjected to an axial force N and a bending moment M, i.e., 

M 

with r=2, simple graphs result for 0 for various values of n and ratios — 
can be obtained. This is accomplished by using 



( 2 . 68 ) 


and finding 


0 = 


r (n) 

e 

sn 


M + dN 

Z 

max 

n 

. M + d*N 

IZ* 

T 

| maxj 



(2.69) 


The starred quantities refer to n=l. These examples are given on Figs. 2-15 
through 2-17. 

To complete the applications of the mathematical theory of the elastic 

analog for non-linear structures some solutions for cases involving rotational 

symmetry (ring, plate, shell of revolution, rotating disk) are tabulated in Tables 

2-1 through 2-4. In these tables a t is the tangential (circumferential) stress, 

a the axial stress and o the radial stress. The same indices are used for 
x , r 

the strain rates e . 

2.7 CREEP BUCKLING 

Design for creep buckling requires that the usual concept of a critical 
(or buckling) load is replaced by the concept of critical lifetimes. The critical 
lifetime is the time necessary for collapse under constant load and tempera- 
ture. The criterion of elastic stability 

H => “ (2.70) 
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Fig. 2-16 -Design Chart for Bending of Rectangular Tube, 
B/H = 0.5, t/H = 0.025 
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Fig. 2-17- Design Chart for Bending of T-Section, b/h = 0.5, t/H 
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Table 2-2 


FORMULAS FOR STRESSES AND CREEP RATES IN PLATES (REF. 2-15) 


Stresses 


Deflection Rate 


Case a 

Simply Support -Distributed 
Load 

W = P 7T R 2 

* n 


l l I X =M 




-2Ro* 




a = - 
r 


ff t = - 


where 


- JyL 2 f(3m+ 1) (l - -4-)] 
\tt mn l V R /J 


3 W r 

2 ~ (3m + 1) - (m + 3) 

8ir mh l a 


. _ 11/1/5 p\-3 l " +, > /s (2ff 4 

W "'" 4 \8 j) n n li- 


i. + 

n+1 


Wf = normal deflection at 
max , 

center 


1 

m = - 


Case b 

Clamped -Distributed Load 
W = on R 2 

~ r\ 


3 ^o- f (3m + 1) \ - (m+1) 
It t mh l R_ 


• -if P.Y 3 ln+l)/t ( 2 « + 
W,„.. - 4 ^24.// ' 


nxi 


cr. = 


3W 


87r mh 


| (m + 3) - (m + 1)| 


Case c 


Simply Supported -Point Load 

J 1 


i, - 7B f ll> V 3< " +n/! < 2 ' 1 +i yn 0 ' 

Wm “ 12 \2tTr ./ ) 
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Table 2-3 

FORMULAS FOR STRESSES AND CREEP RATES IN PRESSURE VESSELS (REFS. 2-15 AND 2-16) 



Stress 

Creep Rates 

Thin -Walled Tubes with 
Internal Pressure 
and Bending 

a t = pr/t 

rr = EI+ MlSL 
CT X 2t 3. 

ITT t 

a =0 
r 

(n-H 

i x = B(ff t - °t + a l > <'VV2 ) 

B(pr/t) n 

Thin-Walled Tubes with 
Internal Pressure and 
Axial Load 

CT t = P r A 

<x x = pr/2t + P/A 

Same as above 

Thin-Walled Tubes with 
Internal Pressure and 
Torsion 

Principal stresses can be determined 
by Mohr’s circle 

±rr-M* 

Principal strain rates can be determined 
and then the strain rates in the axial and 
tangential direction obtained. 

Thick-Walled Cylinder 
with Internal Pressure 

2 2 
R / R v 

rr - n 1 111 ° 

= .^l p ,2.1” / R of 

* p r 2 _ R 2 V 1 

o i 

R i 

a x = P 2 2 

X R - R 

O 1 

2 2 
R / R \ 

° r ‘ P *o - R f ^ ' ' 2 ) 

‘ 4 l(R /B ) 2/n -l “ J ' r ' 

' o' v 

k = 0 

X 

» • 

C = -€ p 

r t 
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Table 2-3 (Continued) 


IV 

I 

00 

oo 



Stress 

Creep Rates 

Thick-Walled Sphere 
with Internal Pressure 
(Ref. 2-15) 

p K 3 / Rq \ 

° t = 3 *3 f 1 + 3 ) 

* R -R7 V 2r / 

0 1 

3 3 

pR , R x 

% ' r 3 -r 3 (*' 7 ") 

O 1 

. B , *n 

e t = 2 < ff t • CT r 

Thin-Walled Cylinders 
Internal Pressure 
(Ref. 2-16) 

°Q = pT ^ 

a x = P r / 2t 

a = 0 
r 

Effective Stress 

S _ SL- 2L 

b _ 2 t 

= (^3/2) n+1 (B) (pr/t) n 
e = 0 

X 

e r = -(^3/2) n ' hl (B) (pr/t) n 



- P r / 2t 

= P r / 2t 

a =0 
r 

Effective Stress 
S = pr/2t 

% = i B <2F> n 

V * i B <tf >" 
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Table 2-3 (Continued) 
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Table 2-4 

FORMULAS FOR STRESSES AND CREEP RATES FOR ROTATING DISKS (REF. 2-15) 
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in which the lateral displacement 6 begins to grow without an increase in load 
P, is augmented by the criterion 


dj5 

dt 


00 


(2.71) 


in which the lateral velocity becomes infinite. 

Numerous creep buckling theories have been devised based on these two 
criteria. They allow the critical lifetime to be computed as a function of load, 
slenderness ratio, initial eccentricity and mechanical properties of the material 
(Refs. 2-17 and 2-18). 


A simple example to illustrate the concept of creep buckling is that of 
the frame shown on Fig. 2-18 in which the concept of a creep hinge described 
by 


9 



(2.72) 


is used. 


P P 



Fig. 2-18 - Frame with Creep Hinge 
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The lateral velocity is 


6 = 0 L 

The bending moment in the creep hinge is 

M = P 6 


(2.73) 


(2.74) 


The differential equation for 6 is given by substituting Eqs. (2.73) and (2.74) 
into Eq. (2.72), 


L 



(2.75) 


With the initial condition 



the solution is given by 

— S' r = (n-l) L (■£—) t 

n -l n-l \U C ) 

o n 

The critical time is reached with 6 becomes infinite 


(2.76) 


(2.77) 


t 

cr 



(2.78) 


Results for the creep buckling of two other structural systems are given in 
Tables 2 -5a and 2 -5b. 


The creep buckling of idealized columns assuming an H-section and 
pin-ended boundaries has been the subject of a number of investigations (Refs. 
2-22 through 2-26 and 2-4). The most popular approach is to assume an initial 
imperfection in sinusoidal form and expand the deformation as Fourier series. 
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Table 2-5a 

CREEP BUCKLING OF BAR WITH "CREEP HINGE" (REF. 2-20) 


p 

1 

creep Hinge 

$1 

1" 


V 

\_i 

Z 

I 

f 



iHlTlAL DEVIATION FR OM STRAf&HT LIME <{> = Zo^ 

DEFLECTION AT t=o <|) C0 = 2. o< ; - Z* a * . _ _p_ 

0 * f 

DEFLECTION AT t 4 1 “ 2<* - Z* 0 6 


D EFi NlTlONS •• 

Euler buckling- 

LOAD 

p , 4 M c 


STRAIN RATE 


♦»= (m/m.) 


ELASTIC STRAIN 


4> e * M/Mc 

Differential 

Equation for 

A (i) 


4<* 

dlt 

\ 1 - — C0S« ] 

L 2 Me j 

f PL N 
\2M h 

j Si'i\ 


CRITICAL Tine WHEW « APPROACHES iNFlVlTr 
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Table 2 -5b 

CREEP BUCKLING OF A SIMPLY SUPPORTED RECTANGULAR PLATE (REF. 2-21) 



real plate 


Equivalent SAubwicH plate 


W 

J.=4._L 

tea 

1_ tT 


I 4 — «— H 


j N IT I AL DEVIATION FROM PLANE 

deflection AT t-o 


el = ~~ for VI a L 


W a W ^ COS Ol COS P 
W(t=o) = w COS o( cos p 


w \ nti 

2w+l ) 


51 <Se 


DEFLECTION AT 

t 


W(t) = 

w (c \t) cos * cos p 

(^E-d 

DEFINITIONS 

EULER 

Buclli'ng- 

STRESS -» 

E firV 

• 


s train 

RATE 

e M = 

kZ* 

-=i W>at«r»^l Constants 


E U LEA 

STRAIN 

e E * 

C?E / E 



EULER 

TIME 

"t E - 

6 e/^« 



D l FFERENTV^L EQUATION FOR VV(t1 Solution ^ 

1. = c, u " C ** ( Cl + ^ C0) 


t E <sL 


l { >* ( C t * w (OX ) 


w li) = 


W ( ° 


critical Time , when w approaches infinity 
— = c, (m f | + C $ 

t E <sL _ ^ L 


VI 

C, 

Ct 

^ 3 

3 

.3 <o<o 

T. b<l 


3*' 

.3 b* 

1.16 

.*120 

S’ 

. n 2 

(.67 

>13 S’ 

7 

. m 

1.36 

.OSU 
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This approach leads to values for t if certain assumptions for n (i.e. f 
n=3, 5, 7, etc. ) are made. No closed form solutions for other than n=odd 
integers have been obtained. 

2.8 NUMERICAL EXAMPLES 

Several numerical examples are presented to demonstrate some of the 
principles discussed in the previous sections. 

Example 1 ; Determine the total strain over a period of time for Rene' 41 
at 1500°F, using the material law of Eq. (2.17). Given are the intercept 
strains as shown on Table 2-6. 


Table 2-6 

INTERCEPT STRAINS (REF. 2-27) 


Stress 

(psi) 

€ 

O 

(in, /in. ) 

35,000 

.00175 

39,500 

.0020 

55,000 

.0030 

55,000 

.0040 

60,000 

.0045 


These strain values include the elastic component. From these data Fig. 2-19 
is constructed. 


The parameters in Eq. (2.17) follow from Eq. (2.18) through (2.20) 


.78 x 10 


-3 


20 x 10 


3 \m 


m 


x 10~ 3 j 

1 l 10 x 

In 2.6 

.955 

fn 2.0 

.693 
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C 


1 



.3 x 10~ 3 

(10,000) 1,38 


.3 x 10~ 3 
3.35 x 10 5 


= 8.95 x 10" 10 


The other two constants are determined from a Lars on -Miller master creep 
curve, where the parameter P according to Eq. (2.2) is 

P = (T + 460) (20 + log t) x 10“ 3 

In Table 2-7 these data for e = 0.2% creep at T = 1500 °F are given for various 
stress levels. 


Table 2-7 

CREEP DATA FROM MASTER CURVE (REF. 2-27) 



Using Eqs. (2.23) through (2.25) the parameters are 


7.4 x 10 


- 5 ' 


2.2 x 10 




/ 40,000 

20,000 


in 


n =• 


i n 33.6 
in 2 


3.515 

.694 


= 5.07 


C„ = 


7.4 x 10 


-5 


7.4 x 10 


-5 


(40,000) 


5.07 


2.24 x 10 


23 


= 3.3 x 10 


-28 
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The equation for the strain at any time t for Rene' 41 at 1500 F is 


o nc in -10 1*38 , , , in -28 5.07, 

e = 8.95 x 10 a + 3.3 x 10 a t 


(2.79) 


Equation (2.79) is evaluated on Fig. 2-20 below for cr=35 ksi. 



Fig. 2-20 - Rene' 41 at T = 1500°F and a = 35 ksi 


Example 2: Determine the creep parameters for Ti-6Af-4V, using Eqs. (2.24) 
and (2.25). Data from Ref. 2-28 show the creep curves given in Fig. 2-21. From 
this figure the following values can be extracted: 


at 


T = 750 F, a = 65,000 psi, 
Ae 


1 


At 


- = ’oX--" = 6.25 x 10" 6 /hr, and 

oU 


at 


T = 750°F, cr = 50,000 psi, 

Ae 2 .05% _ , Q C „ lrt -6 

e 2 = At T T30- - 3 - 85 x 10 /hr • 
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Time (hr) 


Fig. 2-21 - Creep Curves of Annealed MST 6A£-4V 


The parameters of Eq. (2.26) are then 



K 


i_ = 3.85 x 10~ 6 
cr 11 (5.10 4 ) 2 


1.5 x 10 


-15 


To determine (3 in Eq. (2.27) one other data point is needed: 


at T = 800°F (426°C), a = 50,000 psi, 


= (’ 50 = ^TT- = 6. x 10" 5 /hr 


40 
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Then (3 is computed from 



W A x HLi \ 
\ 3.85 X yp_2 
(426°C - 398°C) 


in 15.6 
28 


.0983/°C 


Example 3: Determine the service life of a Ti-6Af-4V bar with axial tension 
load, a = 50 ksi, at elevated temperature, T = 750°F (398°C). 


Cross Section, A 


T 

P 


Fig. 2-22 - Axial Bar 

From Eq. (2.41) it follows for p = 1.0 that the lifetime is 


t 


i 


1 

2 x (5 x 10 4 ) 2 (1.5 x 10" 15 ) 


133,300 hr . 


This service life is large due to the assumption that A/A q = 0 at t = t 
Actually, A/A q « .95 would be a minimum. 
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Example 4: Determine the service life of the bar in Example 3 when a = 50 ksi, 

but T = 750°F +25°F, considering a normal distribution of the temperature 
fluctuations. Then 


AT = 50°F (28°C) 


and the deviation is 



28 

2 


14°C 


and n (t) according to Eq. (2.32) is 




(0.0983) 2 (14) 2 

2 

e 


0.954 

e 


2.6 


The lifetime is 


133,300 

£ ” 2.6 


51,300 hr . 


This is considerably less than at constant temperatures. 


Example 5; Determine the time to creep so that the total tip deflection of a 
tip-loaded cantilever beam is twice the elastic tip deflection. The material 
is Ti-6A£-4V. 



Fig. 2-23 - Cantilever Beam 
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Further data are: T = 750°F (398°C), E T = .76E = .76 (16 x 10 6 ) = 12.16 x 
10 psi, and I = .5 x 1 /l2 = .0417 in. . The elastic tip deflection is from 
Eq. (2.61) 


6 

e 


500 x 10 3 

3 x (12.16 x 10 6 ) (.0417) 


.3285 in. 


From Fig. 2-13 


= _2j2l .5(1) 2 i i _ 

n 2(2) + 1 4 Kz) ’ 

With the creep deflection given, Eq. (2.61) can be solved for the lifetime t^, 


4 i = 


/ PL \ n L Z 

l *n / n+2 


(2.80) 


M K 


m 


For ft = 1.0 (no temperature variation) the lifetime is 


t, = 


(5 X 10 3 \ 2 
\ .0707 / 


.3285 
10 2 


= 1,750 hr 


(1.5 x 


10 “ 15 ) 


For fi = 2.6 (a temperature variation of +_25 F) 


. 1750 Ln-2 

h = TT = 673 hr 


2.9 CREEP ANALYSIS PROGRAM 


The Creep Analysis program is a digital computer program which is 
designed to determine the creep damage for a variable uniaxial loading condi- 
tion. The program will determine the time to rupture, fracture of life expended 
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creep rate, creep strain and total creep strain based on the given material 
data. Three theories of creep deformation are included to predict the total 
creep strain for the load -temperature -time history. These are: 

• Time Hardening Theory (Eq. (2.8)) 

• Strain Hardening Theory (Eq. (2.9)), and 

o Pao-Marin Theory (Eq. (2.12)). 

Plots of creep strain versus time are constructed for these three methods of 
analysis. The program is written in FORTRAN V language and configured for 
execution on the Univac 1108, (see listing in Appendix B). 

The primary difference between the time hardening and the strain harden 
ing theories is the method used for strain accumulation. For the analysis of 
the transient creep range the plots of creep strain versus time are input in 
table form for each stress -temperature level. Interpolation is used to trace 
the strain plot over the time period, At. Both theories sum the change in 
strain values during the time intervals at the various stress levels to arrive 
at the total creep strain value 

n 

e Ae i at At., cr. . 

i=l 

In the change from one stress level to another the time -hardening method 
starts the new time interval at the ending time values of the preceding stress 
level. The strain hardening starts the time interval at the same strain level 
that was reached at the end of the preceding stress level. These two methods 
were illustrated in Fig. 2-11. 

The Pao-Marin theory is an analytical representation of the creep 
strain-time curve by determining parameters based on the material data. 
Parameters are computed for both the transient and minimum rate creep 
strain regions. The strain-hardening method of strain accumulation is used. 
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The other required material input data consist of master creep curves 
such as the Larson-Miller or Manson-Haferd parameters. If both creep rupture 
and percent creep strain curves are input, the service life expended and total 
creep strain are computed. The rupture data are required for the service life 
analysis and the percent creep strain data are required for the creep strain 
calculations. If only one set of data is required or available, the specific 
analysis will be performed and the other bypassed. A maximum allowable 
strain value can also be specified that will terminate the program when the 
total creep strain exceeds this value. This is useful in certain cases where 
excessive deformation may be more critical than rupture time predictions. 

The creep parameter data, both rupture and percent creep, are described 
by inputting specific stress values and their corresponding parameter value in 
table form, starting with the maximum stress value. The program will interpo- 
late between these values in subroutine GIR1 to determine the appropriate param- 
eter for the given stresses. 

As shown in Fig. 2-24 the stress values cr^, etc., and their corre- 

sponding parameter values Pj.P-j.Pj. etc. would be input. For a specified 
stress value f the program interpolates logarithmicly to determine P. Using 
P, the time to rupture or time to creep to a given percent strain can be deter- 
mined. For the creep parameter data the appropriate percent creep strain, 
such as 0.2, 0.5, 1.0% creep, must also be specified for the creep rate computa- 
tions. The Larson-Miller master creep curve for 0.2% creep strain for Rene' 41 
was input in the sample problem in Appendix B. 

The appropriate master parameter curve equation must also be in the 
program for the correct solution. These expressions and their location in the 
program are noted by comment cards, (see program listing in Appendix B). If 
different parameters are used as input data the correct parameter equation 
must be input into the program deck. The original deck is set up with the 
Manson-Succop parameter for the creep rupture data and the Larson-Miller 
parameter for 0.2% creep strain for Rene' 41 data. These parameter expres- 
sions are noted in the program deck. 
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Fig. 2-24 - Master Creep Curve 

A stress level, temperature and time period make up a load case. A 
series of these is constructed to represent the load history of the structure 
to be analyzed. The stress -temperature -time history is idealized as shown in 
Fig. 2-25. 

The time to rupture, fraction of life expended, creep rate and creep strain 
are determined for each of the load conditions. If the stress level falls outside 
the range of stress values given for the parameters a message is printed and the 
program continues to the next load condition. A running sum of the life expended 
and the total creep strain is determined. If the creep strain sum exceeds the 
maximum strain specified, the program terminates. A service life margin is 
computed based on the sum of the fraction parts 


Margin 


1 


1 
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Time 


Fig. 2-25 - Idealized Stress-Temperature-Time History 

where t is the time period and t is the time to rupture at cr. If a negative 
margin occurs the program will terminate. Full output up to the point of 
termination is given if the program is stopped before all load conditions have 
been processed. 

The program solution is carried out in two main DO loops. The first 
loop (DO 131 L=l, 3), specifies the theory to be used. 

L, = 1 Time Hardening Theory 

L = 2 Strain Hardening Theory 
L, = 3 Pao-Marin Theory 

The second primary DO loop (DO 130 1=1, NS), is inside the first loop and con 
trols the calculations for all of the specified NS load -temperature -time condi- 
tions. 
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The program output is given for a sample problem in Appendix B and 
should be self-explanatory with the proper headings. The life expected value 
that is calculated is based on cycling the prescribed load condition and may not 
be appropriate for each problem. The corresponding factor of safety should 
also be used with discretion. It is based on the specified life required and the 
fractional part of the life that is used up and the maximum temperature that 
occurs in the load history. 

The program listing and a sample problem using Rene' 41 data are given 
in Appendix B. Plots of creep strain versus time are generated for each of the 
three strain theories. In the final plot all three theories are compared. The 
time hardening, strain hardening and Pao-Marin plots are denoted by the plotting 
symbols T, S and P, respectively. 
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Section 3 
FATIGUE 

Fatigue analyses must be performed when a structural design includes 
a requirement for a reusable structure with a long lifetime or for a structure 
which must survive a Significant number of fluctuating load cycles. While the 
life characteristics of fatigue sensitive structures should be verified by tests, 
it is also necessary to have the capability to predict the service life with rea- 
sonable accuracy during the design phase. Elevated temperature with the as- 
sociated problems of creep, thermal stresses, temperature and time dependent 
material properties, only complicate the lifetime predictions. 

The elevated temperature fatigue problem can be approached in a num- 
ber of different ways depending on the magnitude of the loads, the temperature 
range, and the number of load cycles. Three of the more important approaches 
to fatigue analysis are: 

• Low-cycle fatigue analysis 

• Fatigue in the creep range 

• Fatigue below the creep range 

Each of these approaches is discussed in the following paragraphs preceded by 
comments on applying factors of safety to fatigue design. 

3.1 FATIGUE - FACTOR OF SAFETY 

It is well known that the fatigue strength of materials exhibits considerable 
scatter and can best be characterized in statistical terms. Unfortunately, few 
materials have been subjected to the vast amount of testing required to quantify 
the statistical properties of the fatigue strength. The engineer is left with a 
need for design criteria which will result in a reliable, fatigue -resistant struc- 
ture without explicitly knowing the statistical characteristics of the material 
properties. 
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Static strength of materials as given in MIL-HDBK-5B have a probabilistic 

basis: 

A basis — The A value is the value above which at least 99% of 
the population is expected to fall, with a confidence of 95%. 

B basis — The B value is the value above which at least 90% of 
the population is expected to fall, with 95% confidence. 

Fatigue data presented in MIL-HDBK-5 are not considered design allowables. 
Instead, fatigue data are presented as typical properties of the material. Two 
kinds of figures are used: (1) individual S-N curve, and (2) constant-life dia- 
grams that can be constructed from a family of S-N curves. These data are, 
therefore, not suitable for estimating survival probability. Other means are 
needed to establish a factor of safety (FS) or probability value to associate with 
lifetime predicted from the curves. 

Inspection of typical S-N curves which show data points or scatter band 

indicate that for a given stress level the spread in the lifetime can be quite 

large, whereas for a given lifetime, the spread in stress is much smaller. 

Figure 1 of Ref. 3-5 shows a set of probability S-N curves for 7075-T6 aluminum. 

This curve is reproduced here as Fig. 3-1. At an alternating stress of 30 ksi 

and a probability of failure of 0.01, the value of N is 1.5 x 10^; for a probability 

8 

of failure of 0.99, N is 1.5 x 10 . The above represents a life factor of 100. 
Between the 0.01 and the 0.50 probability curves there is a life factor of 10. 

These large values of the life factor occur at a low stress level near the en- 
durance limit. At 50 ksi alternating stress, the life factor is only 1.7 between 
the 0.01 and the 0.50 probability curves. 

Inspection of the probability curves reveals that at a specified lifetime, 

*7 

the spread in the fatigue strength is significantly less. At 1.5 x 10 cycles the 

ratio between the 0.50 probability fatigue strength and the 0.01 probability 

5 

strength is 1.25, and at 3 x 10 cycles the ratio is about 1.1. 

Results from fatigue tests of aircraft components are similar to the 
results for fatigue specimens, except that the scatter is even larger for the 
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Fig. 3-1 - S-N Fatigue and Probability Curves for 7075 

Aluminum Alloy (Ref. 3-5) 
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component tests. Figure 3-2 (Ref. 3-7) is based on a large number of com- 
ponent tests from many aircraft programs. Figure 3-2 can be used for 
preliminary estimates of the maximum allowable ultimate design tension stress 
for typical aircraft structures made from aluminum, steel, and titanium alloys. 
The curves indicate the scatter in spectrum fatigue tests and constant amplitude 
fatigue tests of aircraft components. The 50% probability curve is based on 
achieving a fatigue quality comparable to previous aircraft such as the P3V, 
Electra, FI 04, and Model 286 helicopter. The use of design stresses higher 
than the 50% probability design curve must be accompanied by an improvement 
in fatigue quality. Fatigue problems are less likely to occur if design stresses 
below the 50% probability design curve are used. 

The maximum allowable ultimate design tension stress (mean plus alter- 
nating) is expressed as a percentage of the ultimate strength of the material. 
The 50% probability design curve is considered to provide a best estimate of 
the allowable stresses that can be used for the required number of flights or 
cycles. 

Since the S-N curves of Figs. 3-1 and 3-2 are typical of the majority of 
metallic materials, the following conclusions can be drawn: 

1. Fatigue strength data are subject to a large scatter band even 
neglecting complicating factors such as temperature, stress 
concentrations, random loads, etc. 

2. The variability in the fatigue life at a specified stress level is 

much larger than the variability in the fatigue failure stress at 
a given lifetime. ' 

3. The scatter in the predicted service life will increase as the 
applied stress becomes smaller; the scatter becomes extremely 
large as the applied stress approaches the endurance limit. 

4. Additional data are required for most materials before an FS 
can be based on the statistical scatter of the fatigue strength. 
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Fig. 3-2 - Flights or Cycles to Failure 
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3.2 LOW-CYCLE FATIGUE 

When parts are subjected to frequent applications of near-limit loads 
during the service life of a structure, failure can occur due to high stress - 
low cycle fatigue. It is readily apparent from the examination of any typical 
S-N curve or constant-life fatigue diagram that the number of load applica- 
tions which can be sustained is relatively low when high loads are involved. 

High stress-low cycle fatigue is a phenomenon associated with plastic 

strain cycling. While failure may not result in one cycle, sufficient damage 

can accumulate in relatively few cycles to cause rupture. Coffin suggested 

(Ref. 3-10) an empirical relationship between the range of the plastic strain 

(Ae ) and the number of cycles to failure (N). The plastic strain range is 
P 

defined as the difference between the total strain range and the elastic strain 
range. This is illustrated in Fig. 3-3 for a material being cycled. The plastic 
strain range is 



Fig. 3-3 - Cyclic Plastic Strain 
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where 



= Ae 


Act 

E 


(3.1) 


Ae 


is the total strain range and 


Act 

E 


Ae , the elastic strain range. 


The basic relationship for low cycle fatigue failure is 

i 

Ae = CN" 2 (3.2) 

P 

where C is a material constant and N is the number of cycles to failure. For 
many materials the material constant, C, can be found using the elongation, e, 
corresponding to the ultimate tensile stress, F^. Letting N = 1/4 cycle at 
F^ and substituting into the above equation, we have: 

e = C ( 1 /4) ” * (3.3) 


or 


C = e (1/4) 2 = | 

Figure 3-4 shows experimental data (Kennedy Ref. 3-11) for some com- 
mon structural materials which tend to support the empirical relationship. 

Equation (3-2) states that failure occurs when 

N Ae 2 

o- 2 - = 1 (3.4) 

Cf 

Equation (3-4) implies that a cumulative damage law exists, but note 
that it is based on the peak plastic strain in the cycle. This strain can become 
highly localized. Thus, plastic strain computed from gross properties of a 
loaded member could give very unconservative results. This supports the often 
noted situation: Fatigue failures nearly always occur at joints and stress con- 
centrations where the difficulties of analysis are most severe. 
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The discussion above is brief and is only an introduction to high stress - 
low cycle fatigue. The next section, 3.3 Fatigue in the Creep Range, discusses 
low-cycle fatigue in greater detail and relates low-cycle fatigue to thermal 
fatigue and creep failures at elevated temperatures. 

3.3 FATIGUE IN THE CREEP RANGE 

When structures are subjected to combined thermal and mechanical stress 
cycling at elevated temperatures, the mode of fracture can generally be classi- 
fied as. one of two types of failure. If fracture occurs without noticeable de- 
formation it is of the fatigue type, while fracture that is accompanied by notice- 
able deformation is of the creep type. There are several types of combinations 
of mechanical and thermal stresses. In the case where both the thermal and 
mechanical stresses are cyclic and completely reversed the material usually 
fails without deforming significantly. Whereas in a case where the mechanical 
stress is constant with a cyclic thermal input and the material is free to elongate, 
such as a pressurized thick wall cylinder subjected to cyclic temperatures, 
failure would be by excessive deformation. 

In the case of thermal fatigue combined with mechanical mean stress, 
the stress ratio A is used as a parameter. This is the ratio of the amplitude 
of the cyclic thermal stress component to the mechanical mean stress, °’'p/ cr m > 

A simple creep test under steady load and varying temperature corresponds to 
a zero stress ratio, and an infinite stress ratio corresponds to thermal fatigue 
under completely reversed strain cycling. 

The different combinations of thermal and mechanical fatigue will be 
discussed in the following paragraphs. 

3.3.1 Mechanical Fatigue at Elevated Temperatures 

At elevated temperature the fatigue life is greatly affected by the fre- 
quency of stress cycling. This is especially true for high stress amplitudes. 
From the results of Forrest and Tapsell (Ref. 3-15) in Fig. 3-5, the lower 
frequency stress cycling results in shorter cyclic lives to fracture at the 

3-9 

LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D306579 


2 

same stress level. For example, at a stress of 19 to ns /in , the cyclic life 
time at 10 cpm is 12,000 cycles while the corresponding life at 2000 cpm is 
120,000 cycles. 



Endurance-hours to fracture 

A 8000 cpm + 30 cpm 

® 2000 cpm X 10 cpm 

O 125 cpm 


Fig. 3-5 - Fatigue Strength at Various Cycling Speeds 
(From Ref. 3-13) 

Taira (Ref. 3-13) converted the data of Fig. 3-5 into a plot of the plastic 
strain amplitude versus the log of the number of cycles to fracture, N, Fig. 3-6. 



Number of cycles to fracture N 

Fig. 3-6 - Plastic Strain Amplitude vs Number of Cycles to Failure 
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In converting the stress data to plastic strain amplitudes, the dynamic 
stress -strain curves of Fig. 3-7 were used and fracture time was converted to 
number of cycles by taking the frequency into account. 



Alternating strain 

Fig. 3-7 - Dynamic Stress-Strain Curve of 0.17% Carbon Steel at 
500°C (From Ref. 3-13) 

Based on the data of Fig. 3-6, Taira concluded that the fatigue life at high 
temperature is strongly dependent on the magnitude of the amplitude of the 
plastic strain component. Further tests were performed to substantiate this 
conclusion using other materials. In the range of high strain amplitudes the 
plastic strain component of the alternating strain was large and the tests results 
fell on a single curve. In the range of low stress amplitude, where the plastic 
strain component was small, the cyclic life time was found to be more dependent 
on the exposure time, and creep was a factor in the failure. 

To verify the above conclusions, Taira conducted tests at elevated tempera- 
ture for the two materials shown in Figs. 3-8 and 3-9. The figures clearly show 
the divergence in the fatigue life at low strains for the different cyclic rates. 

Based on the fact that the log Ae - log N is a relatively straight line 

P 

(Fig. 3-4 of Section 3.2), a relationship between the increment of damage 

A <j) f caused during one-half cycle of plastic strain Ae was derived 
f P 
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Fig. 3-8 - Strain Amplitude vs Cycles to Fracture for 0.10% Carbon Steel 
at 450°C (From Ref. 3-13) 



Fig. 3-9 - Fatigue Strength of 13 Chromium Steel at 600°C. 

Presented in Strain Amplitude vs Number of 
Cycles to Fracture (From Ref. 3-13) 
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A0 f = MAe p ) u (3.5) 

where \ and u are material constants. It is assumed that failure occurs when 
the sum of the damage increments reaches a critical value 0 • 

0 = 2 N A0 , (3.6) 

o I 

where N is the number of cycles to failure. From the above two relations the 
fundamental equation for fatigue at elevated temperature is 

Ae = (0 /2A)^/ U = constant (3.7) 

p T o' 

Other results have shown that u is very close to a value of 2. When u = 2 
Eq. (3.7) is the same as the equation by Coffin (Ref. 3-10 and Eq. (3-2) of Sec- 
tion 3.2). 

3.3.2 Thermal Fatigue 

In thermal fatigue the strain amplitude Ae and the cyclic temperature 
change between the lower temperature level and the upper temperature level 
are the variables that determine the fracture life. 

The stress -strain hysteresis curve for a temperature cycle is shown in 
Fig. 3-10a and an idealized stress-strain relation for analysis purposes in 
Fig. 3-10b. 

For a simple bar with both ends constrained, the alternating strain in- 
duced by the temperature cycle is 

Ae = Ae + Ae = a (T ? - T.) (3.8) 

0 P Ci x 
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where 


Ae is the elastic strain component 
Ae is the plastic strain component 
a is the coefficient of thermal expansion. 


a 



Fig. 3-10 - Stress-Strain Hysteresis Curve for Thermal Stress 

It is assumed that the material is elastic during the increase in temper- 
ature from T ^ to and also in the cooling from T2 to T^. is a lower 
temperature than and likewise is higher than T^. Plastic deformation 
takes place during the temperature increase from Tj to T2 and during the 
temperature decrease from to T^. 


T 3 = T 2 - e AT 
T 4 = Tj - | AT 


(3.9) 
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where 

| is the ratio of plastic strain component to the total strain. 

Using the concept of cumulative damage similar to that used in the 
previous section and taking X as a function of temperature 

A <f> = X(T)(Ae p ) 2 (3.10) 

Putting this equation in differential form and differentiating gives 


d = 2X (T) 6e • d(6e ) 

Jr 


(3.11) 


Applying the equation to the heating and cooling cycles where plastic deforma- 
tion occurs, Tj - T 2 and - Tj respectively, the total increment of damage 
for one cycle is determined 


A <f> = A <t> + A$ 2 



X (T) (T - T 3 ) dT 


X 

+ 2a^* X (T) (T - T 4 ) dT 


(3.12) 


The coefficient of linear expansion a was assumed constant since the temper- 
ature range is usually small. 


Tests have been performed for repeated thermal cycles varying the 
temperature amplitude AT and also varying the time that the maximum temper- 
ature was held per cycle, Ref. 3-11. The tests show that the value of the max- 
imum temperature reached has more effect on the failure curves than the 
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amplitude of the temperature change. The effect of the time that the maximum 
temperature was applied had considerable effect on the cycle life time as shown 
in Fig. 3-11. Life time is significantly shorter for longer time at maximum 
temperature, leading to the conclusion that failure from thermal fatigue is 
more of a creep rupture than a true fatigue failure. 



Fig. 3-11 - Variation of Life with Time at Maximum Temperature 
(T = 920°C) for Nimonic 90 Under Repeated 

ITlclX 

Thermal Cycles (From Ref. 3-11) 


3.3.3 Thermal Fatigue Combined with Mechanical Stress 


In analyzing the case of thermal fatigue combined with a steady or mean 
mechanical stress the cumulative damage theory used in the previous sections 
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is utilized. It is assumed that failure occurs when the damage accumulated 
by the alternating thermal stress attains a critical value, <f> . In this case 
failure is produced by the combined effect of fatigue and creep damage, <f> £ 
and <f> , respectively. 

<j> = <j> £ + (j> c (3.13) 


Hence failure occurs when <ji = <j> 


K 



i 


From Eqs. (3.6) and (3.10) of the previous sections 


<t> f 2 N X f (T e ) (Ae p ) 2 

6 6 tl_ 

r o r o 


(3.14) 


(3.15) 


where ^f(T e ) is the temperature coefficient of fatigue damage. 


Using the life -fraction theory of cumulative damage for damage in creep, 

a material subjected to a stress cr at a temperature T fails by creep rupture 

at a time t when the damage reaches a critical value, $ . Thus the amount 
r o 

of damage that is absorbed in time, t, is 


t 



o 


The life of a material under creep can be expressed as 

t = X.(T)ct" T 
r c ' ' 


(3.16) 


(3.17) 
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where A (T) and y are constants. Hence Eq. (3-16) becomes 
c 


*c c ° r 

K 7 vn 


dt 


(3.18) 


Substituting Eqs.(3.15) and (3.18) into (3.14) 


N V T e ) * (Ae p )2 , C o 1 

K J M T > 


dt = 1 


(3.19) 


where o is the sum of the mechanical mean stress and thermal stress. 


Using Eq. (3.19) curves similar to those shown in Fig. 3-12 can be con- 
structed. Thermal fatigue tests without mechanical stress and creep rupture 
tests under a constant load and temperature are used to obtain the ordinate 
and absissa, respectively. 



Fig. 3-12 - Stress Fracture Chart for Combined Thermal 
Fatigue and Steady Mechanical Stress (From 
Ref. 3-12) 
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Similarly, design charts can be constructed by using the critical values 
of the alternating stress, o^, that causes fatigue failure, and the mean stress, 
o , that causes creep rupture, each acting separately for the ordinate and 
abcissa respectively. Various combinations of static stress, cr g , and alternating 
stress, o , that will satisfy the interaction relation are assumed (Ref. 3-11). 

cL 


a /<j, + a /o = 1 
a' f s' m 


(3.20) 


Figure 3-13 shows this design chart. Other relations can be constructed 
where the constant p and q depend on the material and tests conditions 


(V 7 / + <W q * 1 < 3 - 21 > 



Fig. 3-13 - Interaction Curve for Combined Thermal 
Fatigue and Steady Mechanical Creep 
Stress 
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Experimental data are needed to help in the construction of such plots 
as those shown in Fig. 3-13. 

3.4 FATIGUE BELOW THE CREEP RANGE 


Fatigue at elevated temperature is dependent on the load intensity, 
temperature, number of load cycles, and the duration of the exposure to 
load and temperature. Below the range where creep becomes a significant 
factor and at moderate loads, standard room temperature design methods 
can be used if the fatigue data are based on appropriate elevated temperature 
fatigue tests. These elevated temperature tests must be performed under 
conditions of temperature, exposure time, and notch sensitivity similar to 
the operating conditions of the full-scale structure. 


While the life characteristics of fatigue — sensitive structures should 
be verified by tests, it is necessary to predict the service life with reasonable 
accuracy during the design phase. The following paragraphs outline methods 
for fatigue design which can be used when creep is not a significant problem. 

The methods for fatigue design which follow are based on procedures 
outlined in Lockheed's Structural Life — Assurance Manual (Ref. 3-7). Methods 
used by the Boeing Company (Ref. 3-8) for fatigue design of the supersonic 
transport (SST) are essentially the same as the Lockheed procedures, with the 
exception of nomenclature. In each case the type of structure or joint is given 
a fatigue rating based on tests and past experience. 


Lockheed uses a fatigue quality index, K, which is somewhat analogous 
to the stress concentration factor, K^.. The fatigue quality index is defined as 
a measure of the ability of a structure to sustain the history of the anticipated 
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service usage. For design the K value of the structure is estimated as dis- 
cussed in Section 3.4.3C. The K value actually achieved in design is based on 
the results of representative fatigue tests. The numerical value of the design 
quality index is defined as the K^. value which yields a D value equal to one, 

(2 n/N = 1.0) for life utilization ratio (D) calculations using the cyclic loads 
sustained in a full scale flight -by -flight fatigue test. The best quality cor- 
responds to the lowest value of K (see Fig. 3-14). 

Boeing uses a detailed fatigue rating (DFR) for components or joints 
which is based on tests and past experience. The DFR is defined as the 
maximum cyclic stress in a constant -amplitude loading cycle at which the 
design detail will withstand 10^ cycles at a stress ratio, R, of 0.06. 

In either case the fatigue quality of the design detail (K or DFR) is 
based on tests and past experience for similar conditions of stress ratio, 
temperature, exposure time, etc. 

3.4.1 General Considerations 

Fatigue analyses are made to provide assurance that the service life of 
a structure will equal or exceed some specified number of service hours 
or flights. To show this, fatigue analyses are made at points in the structure 
using a stress history of anticipated service useage. In addition, fatigue tests 
are usually specified to substantiate the design life requirements. In some 
cases the fatigue test requirements are more stringent than the service life 
requirements. It is the responsibility of the stress engineer to show that 
the structure will meet both the service life and test life requirements. 

When designing a structure to meet the design life requirement, there 
are three main factors to consider. 

1. Choice of material and material processing 

2. Detail design quality built into the structure 
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Failure 


K = 4.5 
Test 



Calculate D values for 
K, =2, 4, 6 and 8 and 

loadings applied in 
flight-by- flight fatigue 
test 


Calculated D Value 


Fig. 3-14 - Illustration of Design Quality Index (K =4.5 Shown). 

Obtained from Flight-by-Flight Fatigue Test Data 
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3. Magnitudes and number of occurrences of constant and alternating 
stresses at critical locations due to the anticipated service load- 
ing conditions. 

The order of these three factors vary, depending on whether the structure is 
basically designed for ultimate strength or fatigue. 

The material and material processing must be chosen with the weight, 
the function of the component and the environmental conditions in mind. The 
choice of material is usually limited because of other considerations, such as 
fracture toughness, corrosion resistance, static strength and stiffness. The 
detail design quality must be kept high by providing gradual changes in load 
paths and by minimizing induced stress gradients. The constant and alternat- 
ing stress levels for normal operating conditions must be kept low and con- 
sistent with the fatigue quality of the structure and the service life requirements. 
Finally high frequencies of fatigue damaging alternating stresses must be 
avoided. The material in the structure must be proportioned to reduce or 
eliminate resonant conditions. 


In selecting a section, on a component or an assembly, for fatigue analy- 
sis, one must consider a critical section, i.e„, a section with high stresses, 
stress concentration, susceptability to fretting, such as threads, lugs, joints, 
etc. However, before conducting a fatigue analysis for a component, the 
following questions should be considered: 

1. Is the component a main load carrying member? If the component 
is a main load carrying member, and is critical in fatigue, then it 
must be representatively fatigue tested. 

2. Is the structural weight of the component a design consideration? 

The material must be selected and the structure designed to meet 
the weight requirements. If the weight of the component is small, 
the component may be purposely designed overstrength and over- 
weight to eliminate potential problems in service. 

3. What is the nature of the normal loading? (Air loads, landing 
loads, centrifugal loads, inertia loads, towing loads, sonic loads, 
temperature loads, etc.?). The loads used in the analysis must be 
as representative as possible of the flight and ground conditions 
anticipated in service. 
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4. Is the component subjected to stress concentrations? Stress 
concentrations occur at notches and section changes along 

the load paths. These stress concentrations must be evaluated 
and taken into account. 

5. Is a limited or unlimited fatigue life required? The methods of 
fatigue analysis and fatigue test substantiation may be different 
for limited life than for unlimited life. 

6. Is the component part of a fail-safe structure? In fail-safe 
structures assurance must be provided by analysis or test 
that damage of a specified size will be detected using normal 
inspection techniques before the damage size becomes critical, 
i.e., catastrophic failure will not occur for fail-safe load condi- 
tions. If the structure is designed to be fail-safe, lower scatter 
factors can be used in the fatigue analysis. 

7. What is the quality of inspection during the manufacture of the 
component ? 

8. What severity of wear and tear is expected in service? If, due 
to wear and tear in service, the fatigue quality or loadings at 

a critical section are affected, this must be considered in the 
fatigue evaluation of the component. 

9. What is the quality of inspection of the component in service? 
The structure at critical sections, like joints, etc., should be 
designed so that the critical sections can be inspected. If this 
is not possible, then a lower design stress must be used. 

10. Is the component repairable? If the component cannot be re- 
paired once a fatigue crack has appeared, it must either be 
designed with sufficient safety margin to prevent cracking or 
it will have to be replaced at regular intervals. 

11. Is the component replaceable? If the component is part of a 
high cost assembly or if the component itself is prohibitively 
expensive to replace, then the component must have sufficient 
safety margin to avoid the necessity for replacement. 

12. Do special environmental conditions exist, such as contact with 
corrosive gases, fluids or solids? The effect of corrosive 
environments must be considered in the fatigue analysis. To 
minimize the detrimental effects of corrosive environments, 
choose materials which are corrosion resistant and/or provide 
for corrosion protection of the material. 

13. Is there a particular type of fatigue problem such as high stress, 
low cycle fatigue or low stress, high cycle fretting fatigue? 


3-24 



LMSC-HREG TR D306579 


14. Will the component be fatigue tested? Fatigue tests are usually 
conducted on major components or assemblies to substantiate 
the fatigue quality of the structure and to establish safe -life re- 
placement times for rotary wing aircraft structure. If fatigue 
tests are not conducted on the components and the structure cannot 
be designed to qualify as fail-safe structure, then lower design 
stresses will have to be used to insure that the structure will have 
an adequate service life. 

3.4.2 Preliminary Selection of Allowable Design Stress 


A. Method of Analysis 


The curve presented in Fig. 3 -15 can be used for preliminary estimates 
of the maximum allowable ultimate design tension stress for typical aircraft 
structure made from aluminum, steel, and titanium alloys. The curves indicate 
the scatter in spectrum fatigue tests and constant amplitude fatigue tests of 
aircraft components. The 50% probability design curve is based on achieving 
a fatigue quality comparable to previous aircraft such as the P3V, Electra, 
F-104 and Model 286 Helicopter. The use of design stresses higher than the 
50% probability design curve must be accompanied by an improvement in fatigue 
quality. Fatigue problems are less likely to occur if design stresses below 
the 50% probability design curve are used. 

The maximum allowable ultimate design tension stress (mean plus alter- 
nating) is expressed as a percentage of the ultimate strength of the material. 
The 50% probability design curve is considered to provide a best estimate of 
the allowable stresses that can be used for the required number of flights or 
cycles. 


To obtain an estimate of the allowable design stresses, the following 
procedure may be used. This procedure should only be used in the early design 
stages of the vehicle before spectra loading data become available . 

1. Determine the number of flights or cycles for which the structure 
must be designed. The number of flights is specified by the design 
specification and/or fatigue and fail-safe policy. 
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Fig. 3-15 - Flights or Cycles to Failure 
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2. Enter Fig. 3-15 with the number of flights or cycles obtained in 
Step 1 and determine the allowable stress as a percentage of 
the ultimate strength of the material. 

3. Multiply the value obtained in Step 2 by the ultimate strength 
of the material. The value obtained is the maximum allowable 
gross area tension stress. 


3.4.3 Determination of Ultimate Design Stresses using an Equivalent 
Ground-Air-Ground Cycle 


During the early phases of design, design ultimate tension stress levels 
must be established for the various structural components. The selection of 
these stress levels must be based on consideration of the anticipated loading 
history to provide assurance that the design life will be achieved. The equiv- 
alent ground -air -ground (GAG) cycle concept provides a relatively rapid means 
for establishing the permissible design ultimate tension stress for preliminary 
analyses of structures subjected to complex spectra of loading. This concept 
may not apply to rotary wing aircraft components where the fatigue loadings 
are not proportional to the ultimate load conditions. 


The equivalent GAG cycle is a measure of the severity of the anticipated 
loading spectra. In concept, it produces the same fatigue damage to the struc- 
ture as would be produced if the complete spectra of loadings were considered; 
and the number of applications of this cycle which can be tolerated by the struc- 
ture is equal to the estimated fatigue life in terms of number of flights or cycles. 
At this time, however, no single definition of the equivalent GAG cycle has been 
shown to be completely adequate. For purposes of preliminary analysis, two 
GAG cycles are considered: 

1. Once per flight peak-to-peak GAG cycle. The loading range for 
this GAG cycle is determined by the maximum and minimum 
loadings which are equaled or exceeded once during each flight. 

2. Average maximum peak-to-peak GAG cycle. The loading range 
for this GAG cycle is the average of the range of the maximum 
and minimum loads which are equaled or exceeded once per flight 
and the range of maximum and minimum loadings which are 
equaled or exceeded once during the life of the structure. 
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The latter value provides a conservative estimate of the equivalent GAG cycle; 
and the two values together envelope the probable range of the cycle. For 
preliminary analyses, the average of the two values obtained is used to deter- 
mine the permissible ultimate design stresses. 

A. Analysis Procedure 

The analysis procedure starts with the spectra of anticipated loadings 
for the structure, the design life, and the required quality index (see Section 
3.4. 3C). 

1. Determine the once -per-flight peak-to-peak GAG loading cycle 
range and the average maximum peak-to-peak GAG loading cycle 
range from the envelope of the anticipated loading spectra. 

2. Calculate the corresponding mean and alternating loadings. 

3. Using the static ultimate design loading (or bending) value and 
a series of selected design ultimate stress values, calculate 

a series of load-to-stress conversion factors. 

4. Calculate the mean and alternating stresses for the two GAG 
cycles for each selected value of design stress. 

5. Enter the appropriate constant-life diagram for the material and 
determine the fatigue life for the two GAG cycles at each value of 
design stress. 

6. Plot the resultant lives as a function of design stress. 

7. Enter the curves at the required design life and obtain the per- 
missible design ultimate stress which is halfway between the 
two curves. 


B. Example Problem 

Determine the permissible ultimate design tension stress for a struc- 
ture subjected to the following conditions: 

Design life: 50,000 flights 
Material: 7075 -T 6 aluminum alloy 
Quality index level: K = 4.0 
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Bending moment spectra envelope: Fig. 3-16 

Ultimate design bending moment; 30 x 10^ in. -lb 


1. From Fig. A-3 

Once -per -flight peak-to-peak GAG bending moment range: 

M = 10 x 10^ in. -lb 
max 

M . = -2 x 10 ^ in. -lb 

min 

Average maximum peak-to-peak GAG bending moment range: 

M = 12.5 x 10 6 in. -lb 
max 

M . = -3 x 10^ in. -lb 

min 


2. Once -per -flight peak-to-peak GAG bending moment cycle: 


M 

= 4.0 x 10 6 in. -lb 

mean 


M ,. = 
alt 

;t 6.0 x 10^ in. -lb 


Average maximum peak-to-peak GAG bending moment cycle: 


M = 4.3 x 10 in. -lb 

mean 

M . = +7.7 x 10^ in. -lb 

alt — 


3. -jj| = .00233 psi/in. -lb @ S deg = 70,000 psi 


= .00200 
= .00167 
= .00133 
= .00100 


= 60,000 
= 50,000 
= 40,000 
= 30,000 
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4. Once-per-flight peak-to-peak GAG stress cycle: 


s m = 934o psi, 

S = 14000 psi @ 

a 

S = 70,000 psi 

0.6 S * 

= 8000 

= 12000 

= 60,000 

= 6670 

= 10000 

= 50,000 

= 5320 

= 8000 

= 4o,ooo 

= 4ooo 

= 6000 

= 30,000 

Average maximum peak-to-peak GAG stress cycle: 


S = 10000 pel, 
m 

S = 17900 psi @ 

£t 

s, = 70,000 psi 

des 

= 8600 

= 15400 

= 60,000 

= 7180 

= 12850 

50,000 

= 5710 

= 10250 

= 40,000 

= 4300 

= 7700 

= 30,000 

5. From Fig. 3 -17 



Once-per-flight peak-to-peak GAG 


N = 9500 cycles 

@ S, = 70,000 psi 


= 22,000 

= 60,000 


= 70,000 

» 50,000 


= 370,000 

= 4o,ooo 


= 3,000,000 

= 30,000 


Average maximum peak- 

to-peak GAG 


N = 5000 cycles 

@ S, = 70,000 psi 


= 8500 

= 60,000 


= 22,000 

^ 50,000 


= 85,000 

= 40,000 


= 600,000 

= 30,000 



\ 
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6. These results are plotted in Fig. 3-18. 

7. From Fig. 3-18, the permissible design ultimate tension stress is 
48,000 psi. 


C. Selection of Fatigue Quality Index, K 

\ 

The fatigue strength of a material is usually evaluated by testing small 
notched or unnotched standard specimens. 


The fatigue strength for a section or a component made from the material 
will normally be lower than the fatigue strength of standard specimens tested 
in the laboratory for a number of reasons. In order to account for the effects 
of conditions not covered by the laboratory tests the fatigue quaity index used 
in the fatigue analysis must include the effects of the following factors: 

Kj. stress concentration factor. This is a factor to be used 
on the nominal local stress at the critical section being 
considered. Its value depends on the geometry at the 
section, i.e., notch effect, change in section on the com- 
ponent, etc. 

K thermal factor allowing for the reduction in allowable 
due to the exposure at elevated temperature. Fatigue 
test data conducted at elevated temperatures are 
usually completed in a relatively short time. If the 
material will be exposed to elevated temperature for 
extensive periods of time in service, then the strength 
degradation effects due to temperature exposure must 
be accounted for. In lieu of data for the correct ex- 
posure conditions the K x correction factor can be used 
where K x is equal to the ultimate strength of the mate- 
rial at temperature for a short period of time (usually 
2 hour exposure) divided by the ultimate strength of 
the material after exposure to temperature for the 
correct period of time. In order to substantiate the 
values used, fatigue tests should be conducted under 
the correct exposure conditions. 

K corrosion fatigue factor accounting for the effect of 
corrosion under service conditions. In some cases 
the environmental conditions may be beneficial, as for 
example vacuum as compared to air, whereas the more 
usual environmental conditions will be detrimental as 
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compared to tests conducted in air. Corrosion factors 
are generally not available. Because the fatigue strength 
can be greatly reduced by corrosion, corrosion must be 
avoided whenever possible. 

Other effects, such as fretting, joint eccentricity, differences in heat 
treatment, etc., should also be accounted for. 

The above factors should be multiplied together so that the fatigue quality 
index, K, to be used in conjunction with the nominal local fatigue stress in the 
fatigue analysis becomes: 

K = K x K x K (3.22) 

t x c 

Reduced S-N data can be used in lieu of a K factor if the percent reduction is 
based on representative fatigue test data. 

The design goal is to achieve the lowest practical fatigue quality index. 
However, the minimum K values that can be used for fatigue analysis at Lock- 
heed Aircraft Corporation are specified in the table below. 


Minimum K Values for Fatigue Analysis 

Structures that will be repaired if damaged in 
service such as shell structure 

Structures that will be replaced rather than 
repaired in service such as landing gear 
structure, dynamic rotating components, 
control system components, etc. 


4.0 Aluminum 
4.5 Titanium 


2.0 


The minimum K value of 4 is based on analysis of previous service experience 
and test data. Early service failures are characterized by K values greater 
than 4 whereas parts which have demonstrated adequate service life invariably 
give K values less than 4. In this manner successful past service history is 
projected into new design. This comparative base is most important to 
maintain in new design. 
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K values less than those given in the preceding table should be used 
cautiously. For fatigue critical structure, fatigue quality must be sub- 
stantiated by test of the structure. 

3.4.4 Palmgren-Miner Method of Fatigue Analysis 


A. Method of Analysis 


The method of analysis presented in this section is based on the Palmgren- 
Miner theory of linear cumulative fatigue damage. While it is recognized that 
this method of analysis is not precise, it will yield reasonable estimates of the 
service life of structural components when used as described in this section. 

The basic equation is expressed as follows: 


D 


n, 


N. 


n„ 


N. 


+ . 






(3.23) 


where D = calculated life utilization ratio 

th. 

n. = number of loading cycles applied at the i stress level 

= number of loading cycles to failure for the i n stress level from 
the relevant constant-life diagram. The relevant constant -life 
diagram is the one which applies to the material and fatigue 
quality index of the section under consideration. 


n. 

= cycle ratio 


k 


number of stress levels considered. 


Use AS = 3 to 4% of F, . 
a tu 


The method of analysis using Eq.(3.23) above as a basis, follows: 


1. Select a fatigue quality index K, for the section under considera- 
tion, as described in Section 3.4. 3C. 
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2. Obtain a loading spectra for the section under consideration. 

These spectra should cover all the ground cases and all the 
flight cases including the effect of maneuver, gust, etc. A 
constant and an alternating load should be given for each case, 
together with the anticipated number of loading cycles for the 
design life of the vehicle. 

If the loads are given in graphical form as cumulative loading 
spectra, then discrete loading distributions must be developed. 

3. Convert the loads, for each case, to stresses at the section 
under consideration. 

4. Obtain a constant -life diagram for = K, for the material 
under consideration. 

5. Calculate the cycle ratio for each loading case, including the 
effect of the GAG cycle, and add the cycle ratios to obtain 
the life utilization ratio "D." 

6. Calculate the life in hours for the section under consideration 

from R,L, 

L c = -V 1 < 3 ' 24 > 

where L^ = the life span in hours represented in the spectra 
used in the analysis. 

R2 = a reduction coefficient used to assure a specified prob- 
ability of obtaining a test life equal to or greater than the 
calculated life. Definitions of the GAG cycles used for fatigue 
analysis are illustrated in Fig.3-19V Values of R£ for various 
probabilities are given in the table below. The values in the 
table are those used by the Lockheed -California Company for 
design of aircraft structures. The values of R 2 are based on 
experience with typical aircraft structures and are presented 
for information and not as firm recommendations. 


Test Life Reduction Coefficient, R£ 


Probability of Conservative 
Prediction of Test Results 


For (GAG)pp For Spectrum 

of GAG 


50% 

90% 

95% 


.50 

.60 

.25 

.33 

.20 

.25 
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Stress 



Fig>3-19 - Definition of GAG Cycles Used in Fatigue Analysis 
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7. When the value of the fatigue quality index n K !l is not known, 
evaluate the life utilization ratio and hence the calculated 
life for several values of n K, M as described in Steps 1 to 6. 

Plot K versus the calculated life as shown in the plot below. 


:V 




1 L.J, LLLLL 

1 L-L. LJLLJi 

' 

i . LI.I1H1 

l I i i n it 


°icf 1C) 4 10 5 10 c 10 7 

Calculated Life — Hours 

Fig. 3-2 0 - Relation Between Fatigue Quality Index and Calculated Life 


The calculated life can thus be obtained for any value of M K n 
within the range considered. Hence, if the range has been 
chosen correctly, the calculated life for the section under 
consideration can be obtained at a glance once !, K n has been 
established, as discussed in Section 3.4.3C. 

8. If the calculated life is less than the design life, the following 
steps should be taken in order of priority: 

a. reanalyze the part using more stress levels to represent 
the loading spectra for cases which have large cycle ratios. 

b. redesign the part to bring the calculated life up to or above 
the design life. 

c. conduct fatigue tests to demonstrate that the component 
will meet the design life requirement. The minimum test 
life required to demonstrate the design life is a function 
of the type of test, number of tests, type material, and 
magnitude of the stress levels. 
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9. Compliance with the design life goal must be demonstrated by either 

a. Representative fatigue tests of the actual part or of parts with 
similar materials and stress concentrations. 

or b. Previous service experience of parts of similar structural 
design taking into account the differences in structural 
design and operating conditions and procedures. 
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Section 4 

PROCEDURES FOR STRUCTURAL ANALYSIS AT 
ELEVATED TEMPERATURE 

Of primary importance and the first step in any structural analysis at 
elevated temperature is the accurate prediction of the loads and temperature 
environment. Any conservatism in the predicted loads and temperatures re- 
sults in an unwarranted weight penalty in the structural design. Accurate 
temperature predictions are extremely important as the temperature approaches 
the upper temperature limit of the material. At the upper limit the strength 
properties decrease sharply, the thermal stresses increase and creep be- 
comes a major consideration in the design. 

Most elevated temperature structural design problems are too complex 
to simplify to handbook or hand calculations without undue conservatism being 
added to the analysis. Finite element or other mathematical computer models 
are usually constructed to assist in the analysis. These models offer the 
capability of varying the material properties, such as E and a, from element 
to element to correspond to the thermal distributions. This allows a more 
accurate simulation of the real structure as greater detail can be provided 
in problem areas. This method is accurate for thin-skin and/or structures 
approaching the plane-stress state. Based on the results of the initial com- 
puter run the math model can be further refined by inputting reduced modulus 
values in the elements that have been shown to exceed the proportional limit. 
This procedure has been used to analyze complex structures with large ther- 
mal gradients, Refs. 4-1 and 4-2. 

In determining the failure stress for the structural components, based 
on what is considered failure (buckling, plastic yielding, excessive deforma- 
tion, etc. ),the complete temperature and load history must be covered. A 
flange that might fail by buckling at a room temperature environment could 
fail by plastic yielding at an elevated temperature environment. A good 
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practice to follow is to construct a plot of critical stress versus temperature, 
Fig. 4-1, so that the minimum value will always be used at the appropriate 
temperature. 



Fig. 4-1 - Critical Stress 

4.1 MATERIAL ALLOWABLES AT ELEVATED TEMPERATURE 

The material properties used in the analysis for most metallic aerospac 
structures are taken from MIL-HDBK-5 (Ref. 4-3). The room temperature de 
sign allowables are presented on one of four bases: 

A basis - The value above which 99% of the population of values is ex- 
pected to fall with a confidence of 95%. 

B basis - The value above which 90% of the population of values is ex- 
pected to fall with a confidence of 95%. 

S basis - The S value is the minimum value specified by the governing 
group responsible for establishing material specifications. The statistical 
assurance associated with this value is not known. 
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Typical basis - The typical value is an average value with no statistical 
assurance associated with it. 

Due to the amount of data required, usually only tensile ultimate and 
yield strengths are determined on A or B basis. Also, only tensile ultimate 
and yield strengths and elongation are specified and termed S values. Ratioing 
procedures have been established with which other property values are computed 
to have approximately the same assurance levels as the values for tensile 
ultimate and yield. All elastic modulus values, Poisson's ratio values, and 
physical properties are presented as average values. A basis values are always 
used unless the specified value, S value, is lower. Appropriate footnotes are 
given in MIL-HDBK-5 for these cases. 

Elevated temperature properties are usually presented graphically as a 
percentage of the room temperature value and the effects of time at temperature 
are included. If the room temperature value has an A basis the elevated tempera 
ture values are assumed to have an A basis also. 

The stress-strain and stress -tangent modulus curves in MIL-HDBK-5 
are "typical" curves. Typical curves means that the stress -strain data have 
been adjusted to reflect average values of the elastic modulus and typical values 
of the 0.2% offset yield strength in tension and compression. Where stress- 
strain curves at elevated temperature levels are given and A basis values are 
needed, the following method can be used to obtain representative values. It 
is necessary to reduce the stress values from the typical stress strain curve 
to make them compatible with the A base table values. The F. value given 
in the table is ratioed to the 0.2% offset stress from the curve to obtain an 
appropriate working value. 


R 


= V f 0. 


2 % 


(4.1) 


f = R • a 


(4.2) 
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Any stress value (cr) read from the stress-strain curve above the proportional 
limit is factored by this ratio to make it compatible with the MIL-HDBK-5 
tabular data. 



Strain 

Fig. 4-2 - Elevated Temperature Allowable From Typical Stress- 
Strain Curves 

4.2 STRESSES IN THE INELASTIC RANGE 

For the general case of an indeterminate structure subjected to mechanical 
loads, a nonlinear analysis is required for accurate computation of the load- 
deflection response. The statically indeterminate structure subjected to 
mechanical loads differs from both the statically determinate structure and 
the indeterminate structure with thermal loads: 

1. The internal loads in the statically determinate structure are 

independent of the deformations unless the deformations increase 
to the point where geometric nonlinearities become significant. 
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2. The internal loads of an indeterminate structure subjected to 
thermal loading tend to be self-limiting as the structure goes 
into the inelastic range, i.e., a small amount of plastic de- 
formation relieves the thermal stresses, 

3. Inelastic deformation of an indeterminate structure subjected 
to mechanical loads only redistributes the internal loads in 
the structure. 


Since most finite element and other general structural analysis programs 
are linear, the stresses and strains in the elements are always computed using 
the linear relationship 


e = f/E (4.3) 

For statically determinate structures in the inelastic material range, the 
stresses computed by a linear analysis are correct and the strain can be de- 
termined by referring to the stress -strain curve. 

When a structure has only thermal loads and is operating in the nonlinear 
range of the material, the thermal stresses can be approximated by assuming 
the strain computed by a linear analysis is correct. The thermal stresses are 
determined by referring to the stress-strain curve or by using an analytical 
expression for the stress -strain relationship such as the Ramberg -Osgood 
model. This approximation probably will be sufficiently accurate if the material 
is not operating too far into the nonlinear range. The implicit assumption of 
this approach is that the deformed shape of the structure can be approximated 
by the deformed shape of a linear structure. For structures operating well into 
the nonlinear range a nonlinear analysis must be used as for any statically 
indeterminate structure. 


4.3 RECOMMENDATIONS 


1. Use a temperature margin AT that is added to the limit temperature to 
account for temperature prediction unknowns and temperature overshoot 
in variable temperature environments. This temperature margin could 
be a certain temperature value, AT = T, or a percentage of the tempera- 
ture variation in the structure, AT = k (T - T 

max amb' 
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2. Development of non-linear structural analysis programs that will take 
variation of material properties as a function of temperature and time. 

3. Further material tests at elevated temperature levels to determine 
material properties on an "A" basis instead of the present typical 
properties. 
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Section 5 

DESIGN METHODS FOR COMBINED ENVIRONMENTS 

In the previous sections methods for structural design at elevated tem- 
perature have been discussed considering creep, creep rupture, fatigue, 
thermal stresses, and methods of stress analysis. This section discusses 
methods of combining load -temperature -time environments for structures 
operating at elevated temperatures. 

To define quantitatively the confidence in structural integrity or the 
reliability of the structural design, one must rely on probability theory. 
Unfortunately the statistical approach is cumbersome for use during the design 
phase, and statistical data are seldom available when the structure is being 
designed (often adequate statistical data never becomes available). 

The alternative to the pure statistical approach to structural reliability 
is to apply design factors, similar to room temperature factors of safety, to 
the mechanical loads, thermal loads or temperatures, and the service life. 

The structural reliability cannot be quantitatively determined using this design 
factor approach; but, if the design factors are applied consistantly with sound 
reliability principles, confidence in the structural integrity can be established. 

Although the general principles of elevated temperature design are appli- 
cable to any type structure, the following sections are slanted more to a flight 
vehicle such as the Space Shuttle. 

5.1 DESIGN CONSIDERATIONS 

A primary design consideration in flight structures is minimum weight 
in addition to costs and manufacturing difficulties. When the materials are 
selected for high temperature application, curves of efficiency parameters 
for modulus, strength and stiffness (stability) are studied first (Ref. 5-1). For 
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a given configuration, the sensitivity level to high temperatures must then be 
established. This should give an estimate of: (1) special elevated temperatures 
design and testing procedures that may be required; (2) if the usual room- 
temperature design approach is acceptable; or (3) if the configuration is in- 
admissible. 

An acceptable design will avoid both structural and functional failure. 
Structural failure is characterized by a disintegration of a part of the structure 
(fracture, rupture). Functional failure is the result of large permanent defor- 
mations (plastic or creep strains). In considering creep, a design goal is to 
make the creep rate extremely low and the time to fracture extremely long. 

The concept of damage accumulation is essential in design for service life. 
Structural damage under low intensity loads of long duration reduces the re- 
sistance to loads of high intensity and short duration. 

Therefore, the designer must consider the following types of structural 
failure: 


1. Failure of the undamaged structure due to short duration 
loads having an extremely small probability of occurrence 
(ultimate load failure). 

2. Failure of the damaged structure by ultimate loads of a 
somewhat smaller amplitude but higher probability of 
occurrence, 

3. Functional failure due to excessive deformations from 
creep or plastic strain. 

In designing the structure to resist the above failures it is necessary to 
consider the interaction of the various loads (sustained and alternating, thermal 
stresses, etc.) and the temperature environment as functions of time. Numerous 
methods of estimating the structural response to the combined environments 
have been proposed (Refs. 5-2 through 5-6), but none of these methods are 
wholly satisfactory for all cases of loads, temperatures, and materials. In 
most cases elevated temperature testing which simulates the material load- 
temperature-time environment is required. 
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5.2 INTERACTION CURVES 

In Section 3 the interaction between thermal fatigue and mechanical 
stress was discussed. Similar interactions exist between mechanical fatigue 
and creep stresses at elevated temperature. 

Because damage does not accumulate at very low values of sustained stress 
the effect of creep damage can be neglected for small values of the ratio of sus- 
tained to alternating stress. Similarly, very small alternating stresses have 
little effect on the creep response of elevated temperature structures. For 
these reasons elliptic interaction curves are sometimes proposed as a first 
approximation (Ref. 5-2). Examples of the elliptic interaction curves are 
shown in Fig. 5-1. These curves show a fracture or damage surface as a 
function of temperature, T, and a time, t. The creep time, t, and fatigue life, 

N, are related by a cyclic frequency, u: 

t = N/w 

Note that these curves are not the same as the modified Goodman diagrams at 
normal temperature which relate the mean and alternating components of the 
fatigue life. The creep-fatigue interaction curves can only be used at elevated 
temperature when creep damage accumulates under the sustained stress. If 
a number of these curves are available at different lifetimes, t, a three- 
dimensional plot can be constructed as shown in Fig. 5-2. In the figure, the 
intersection of the failure surface with the c^-t plane is the S-N curve at tem- 
perature T. The intersection with the a -t plane represents the creep rupture 
curve. 

Figures 5-3 through 5-9 are from MIL-HDBK-5B and represent material 
test data for creep-fatigue interaction. In Figs. 5-3 and 5-4 the constant life 
diagrams are typical of creep-fatigue data from MIL-HDBK-5B (Ref. 5-8). 

Figure 5-3 presents fatigue-creep rupture and fatigue-0. 2% creep data, while 
Fig. 5-4 shows fatigue -creep rupture and fatigue -0.5% creep data; both figures 
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Fig. 5-1 - Elliptic Interaction Curves for Failure 
at Various Temperatures 



Fig. 5-2 - Interaction Surfaces for Failure 
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A* 4.0 2.33 1.5 I 0.67 0.43 0.25 0.11 

R * -0.6 -04 -02 0 0.2 04 06 08 I 



Minimum Stress, ksi 


FIGURE 6.3.7.2.8(b). Typical constant-life diagram for fatigue and dynamic creep behavior of solution treated and aged 
M-252 forgings at 1500 F. 


Correlative Information for Figure 63.7.2.8(b) 

Product Form : Forged bar, 1-3/4-inches diameter 

Properties: TUS, ksi TYS.ksi Temp,F 

176.0 100.0 RT 

100.0 73.0 1500 F 

Specimen Details : Unnotched 

0.250-inch diameter 

Surface Condition: Longitudinally polished with 240, 400 and 600 grit belts to provide surface finish of 5-8 
RMS. 

Heat treatment included solution treatment at 1950 F for 4 hours, air cooled: aging at 
1400 F for 1 5 hours (packed in cast iron chips), air cooled. 


Fig. 5-3 - Typical MIL-HDBK-5B Creep-Fatigue Data 


Test Parameters : 

Loading - Axial 
Frequency - 1800 cpm 
Temperature - 1 500 F 
Atmosphere - Air 
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FIGURE 6.3.7.2.8(a). Typical constant-life diagram for fatigue and dynamic creep behavior of solution treated and agedM-252 
forgings at 1500 F. 


Correlative Information for Figure 63.7.2.8(a) 

Product Form: Forged bar, 1-3/4-inches diameter 

Properties: TUS, ksi TYS, ksi Temp, F 

176.0 100.0 RT 

100.0 73.0 1500 F 

Specimen Details: Unnotched 

0.250-inch diameter 

Longitudinally polished with 240, 400 and 600 grit belts to provide surface finish of 5-8 
RMS. 

Heat treatment included solution treatment at 1950 F for 4 hours, air cooled: aging at 
1400 F for 15 hours (Packed in castiron chips), air cooled. 


Fig. 5-4 - Typical MIL-HDBK-5B Creep-Fatigue Data 


Surface Condition: 


Test Parameters: 

Loading • Axial 
Frequency - 1800 cpm 
Temperature - 1 500 F 
Atmosphere - Air 
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being for M252 forgings at 1500F. The quasi-elliptical shape of the interaction 
curves are apparent. 

Figures 5-5 and 5-6 present fracture interaction data for Udimet 500 
alloy at 1200F and 1650F, respectively; the solid lines are for unnotched 
specimens and the dashed lines for notched speciments. Note that, while the 
notches greatly reduce the fatigue strength, creep rupture properties are im- 
proved for a sustained stress acting alone. 

Figures 5-7 through 5-9 are from MIL-HDBK-5B for N-155 bar stock. 
These figures are similar to Fig. 5-3 through 5-6 except that the data are 
presented on a cycle-to-failure basis instead of a time basis. As the frequency 
of loading is given, the conversion to a time scale for creep is apparent. 

Material data of the form shown in Figs. 5-3 through 5-9 can be used 
to generate interaction surfaces of the form shown schematically in Fig. 5-2. 
These interaction curves will be valid for a particular loading spectra only 
if the frequency term id is correct for relating creep time and fatigue cycles. 

In the absence of creep-fatigue interaction data, an approach similar to 
that discussed in Section 3.3.3 and equations (3-13) through (3-21) can be used. 

5.3 DESIGN FACTORS AT ELEVATED TEMPERATURE 

As stated previously the goal of the designer is to arrive at a structure 
which has a high probability of survival (reliability) without excessive penalties 
in terms of weight and cost. A probabilistic approach to structural reliability 
is desirable if adequate information is available; unfortunately, this is seldom 
the case during the design phase of new hardware. Since it is not presently 
practical to use a statistical reliability approach to the design of structure, 
an alternative procedure is to use reliability principles to arrive at a rational 
approach for applying design factors. 
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Test Conditions 


Unnotched 


Notched 

K t = 3i4 

A*a> A* 0.25 


100 IZO 140 160 180 200 


Minimum Stress, Ksi 


FIGURE 6.3.9.1 .8(a). Typical constant-life diagram for fatigue behavior of solution-treated and aged Udimet 500 alloy 
bar at 1200 F. 


Correlative Information for Figure 63.9.1 .8(a) 

Product Form: Rolled bar, 3/4-inch diameter 

Properties : TUS, ksi TYS, ksi VTemp, F 

(no properties Given) 

Speciment Details: Unnotched Notched, V-Groove,Kt -=3.4 

0.200-inch diameter 0.375-inch, gross diameter 

0.250-inch, net diameter 
0.010-inch, root radius, r 
60° flank angle, co 


Test Parameters : 

Loading - Axial 
Frequency - 3600 cpm 
Temperature - 1200 F 
Atmosphere - Air 


Kfl = 2.41, p = 0.0022 inch, where = 1 + 


Surface Condition : Unnotched : Longitudinal polish with 400 grit. 

Notched: Notched polish with 600 grit. 


K t -1 


it p 


7T-00 M r 


Fig. 5-5 - Typical MIL-HDBK-5B Creep-Fatigue Data 
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A= 4.0 
R= -06 

2.33 

-0.4 

1.5 

-0.2 

1 

0 

0.67 

0.2 

0.43 0.25 0.11 0 

0.4 06 08 IjO 



Minimum Stress, ksi 




FIGURE 6. 3. 9.1 .8(b). Typical constant-life diagram for fatigue behavior of solution-treated and aged Udimet 500 alloy 
bar at 1650 F. 


Correlative Information for Figure 6.3.9.1.8(b) 


Product Form : Rolled bar , 3 /4-inch diameter 

Test Parameters: 

Properties: TUS, ksi TYS, ksi Temp F 

84.0 - 1650 F (Unnotched) 

1650 F (Notched) 

Loading - Axial 
Frequency - 3600 cpm 
Temperature - 1650 F 
Atmosphere * Air 


Specimen Details: 

Un notched 
0.200-inch diameter 

Notched, V-Groove, K t = 3.4 

0.375-inch, gross diameter 
0.250-inch, net diameter 
0.010-inch, root radius, r 
60° flank angle, ou 

Surface Condition: 

Unnotched: Longitudinal polish with 400 grit. 
Notched: Notch polish with 600 grit. 


Fig. 5-6 - Typical MIL-HDBK-5B Creep-Fatigue Data 
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FIGURE 6.2.2.1.8(a). Typical constant-life fatigue diagram for solution-treated and aged N-155 bar at 1200 F. 


Correlative Information for Figure 6.2.2.1.8(a) 


Product Form : Rolled bar, 1-inch diameter 

Properties: TUS,ksi TYS,ksi Temp,F 

80.0 - 1200 F 


Specimen Details: Unnotched 

(L225-inch 


Notched, V-Groove, K t = 5.1 
0.319-inch, gross diameter 
0.22 5 -inch, net diameter 
0.005-inch, root radius, r 
60° flank angle, co 


Test Parameters : 

Loading - Axial 
Frequency - 1500 cpm 
Temperature - 1200 F 
Atmosphere - Air 


Kj 4 = 3.37, p - 0.0012 inch, where = 1 + 



Surface Condition : Unnotched specimens were longitudinally polished with 400 grit paper. Notched specimens 
were lathe turned in the notch with a carbide tool. 

Heat treatment involved solution treatment at 2200 F for 1 hour, water quench; aging 
treatment at 1400 F for 16 hours. 


Fig. 5-7 - Typical MIL-HDBK-5B Creep- Fatigue Data 
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FIGURE 6.2.2. 1 .8(b). Typical constant-life fatigue diagram for solution-treated and aged N-155 bar at 1200 F 


Correlative Information for Figure 6.2.2.1.8(b) 


Product Form: Rolled bar, 1-inch diameter Test Parameters: 


Properties: TUS, ksi TYS, ksi 

Temp, F 

Loading - Axial 
Frequency - 1500 cpm 

80.0 

1200 F 

Temperature - 1200 F 

Specimen Details: Unnotched 

0.255-inch diameter 

Notched, V-Groove, K t = 1 .73 
0.319-inch, gross diameter 

Atmosphere - Air 


0.225 -inch, net diameter 
0.050-inch, root radius, r 
60° flank angle, co 

Notched, V-Groove, K t = 3.2 

0.010-inch root radius, r 

other dimensions are as given above 



Surface Condition: Unnotched specimens were longitudinally polished with 400 grit paper. Notched specimens 
were lathe turned in the notch with a carbide tool. 

Heat treatment involved solution treatment at 2200 F for 1 hour, water quench; aging 
treatment at 1400 F for 16 hours. 


Fig. 5-8 - Typical MIL-HDBK-5B Creep-Fatigue Data 
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Minimum Stress, ksi 

FIGURE 6.2.2.1.8(c). Typical constant-life fatigue diagram for solution-treated and aged N-l 55 bar at 1500 F 


Correlative Information for Figure 6.2.2.1.8(c) 
Product Form: Rolled bar, 1-inch diameter 


Properties: TUS.ksi 

TYS.ksi 

Temp, F 

50.0 

- 

1500 F (Unnotched) 



1500 F (Notched) 


Specimen Details: Unnotched Notched, V-Groove, K t = 5.1 

0.225-inch diameter 0.319-inch, gross diameter 

0.225-inch , net diameter 
0.005 -inch, root radius, r 
60° flank angle, co 


Test Parameters: 

Loading - Axial 
Frequency - 1500 cpm 
Temperature - 1500 F 
Atmosphere - Air 


Kj*j = 2.16, p = .016 inch, where = 1 


+ 


K t -1 



Surface Condition: Unnotched specimens were longitudinally polished with 400 grit paper. Notched specimens 
were lathe turned in the notch with a carbide tool. 

Heat treatment involved solution treatment at 2200 F for 1 hour, water quench; aging 
treatment at 1400 F for 16 hours. 


Fig. 5-9 - Typical MIL-HDBK-5B Creep-Fatigue Data 
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This section discus ses the rationale of applying design factors to a 
structure operating at elevated temperatures. For purposes of the discussion 
it is assumed that failure surfaces in terms of load, temperature, and time 
can be defined; and that reliability values can be associated with the applied 
loads and temperatures. 

Consider Fig. 5- 10a which shows a failure surface for a structural com- 
ponent in terms of the normalized failure load, P^/p^i the normalized lifetime, 
t/t ; and the component temperature, T. Also shown on Fig. 5-10a are the 
normalized design load, P-q/p^J the expected or limit temperature, T^; and the 
maximum or ultimate temperature, T . The loads are considered to be the 
loads applied to the structure, and the loads and stresses in the component 
are proportional to the applied loads only for a linear, elastic structure. 

p^ is the allowable or failure load as a function of (T, t) 
p^ is the limit or maximum expected load in service 
P D is the design load and 

p D = K • p* 


where 

K is a design factor of safety on the limit load. 

T^ is the limit or maximum expected temperature on the 
component. T^ is the temperature in the component 

associated with the limit conditions of heating, Q^(t). 

T^ = T^ + AT, the ultimate temperature condition in the 

component where the AT term accounts for unknowns 
in the heating environment and errors in computing 
the component temperature. 

t is the service life of the structure and is associated 
with the number of flights or cycles by a loading 
frequency, a). 
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5.3.1 Reliability Considerations for Design Loads and Temperatures 

The undamaged structure should resist ultimate load failure due to a 
load of a very low probability of occurrence, and the damaged structure should 
resist failure due to ultimate loads of a somewhat lower probability of occur- 
rence. Structural reliability should be comparable to that given by the factor 
of safety on mechanical loads at room temperature. 

In room temperature design the limit loads are computed with some 
(generally unknown) probability that the loads will not be exceeded in normal 
service. The design load is determined using a factor of safety, which, based 
on experience, gives a reliable structure. The cumulative probabilities for the 
loads are 


p (p^) = probability (p < p^) 
p (P u ) = probability (p < K Q • p^ = pj 


where K.^ is the design factor of safety. 

The complement to the cumulative probability is the probability of 
exceedance: 

R (P^) = 1 " p (Pj 0 ) = probability (P > P^) 


and 

P(p u ) > P( Pj! ) 

R (p u )< R (P i ) 

The product rule of probability theory states that the probability that two inde- 
pendent random events will occur concurrently is the product of the individual 
probabilities of occurrence. If the mechanical load and the temperatures are 
independent random events the joint probability of exceedance can be expressed 
as 

R(p 1 ,T 1 ) = R(pj) • R(Tj) 
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which is the joint probability that: 

(P >pj) and (T > Tj) 

Consider Fig. 5- 10b. Loads analyses predict a maximum limit load that 

the structure will experience and this load has an unknown probability, P(py. 

Experience has shown that for room temperature analysis acceptable reliability 

results when a specified factor of safety is applied to the limit load to obtain a 

design load. The design load has an unknown but acceptable probability of 

P(Pj- ) ). Similarly, probabilities can be associated with the temperatures, 

and T . 
u 

The joint probabilities of exceeding a stress and temperature can be 
estimated if the probabilities of the stresses and temperatures occurring are 
written 

P( Pi ) = i - R(py 

P(p d ) = 1 - R(P d ) 

pay = i - Ray 

P(T u ) = 1 - R(T U ) 

and 

R(p d ) < R(Pj) 

R(T U ) < Ray 

Consider the joint probability of exceeding the four numbered points on Fig. 5- 10b. 

R i = R(p i* T f ) = R(py*Ray 
R 2 = R(p d ,t y = R(p D ). Ray 

R 3 = R (P D ' T u ) - r «Pd> •*<*»> 

R 4 = R(p r T u ) = R(py.R(T u ) 
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Table 5-1 is an illustration of the joint probabilities of exceeding the 
stress and temperature for the following condition: 


R( Pjf ) = Rf Tj) = 0.10 
R(p d ) = R(T U ) = 0.01 


P( Pi ) = P(Tj) 
P(P D ) = P(T 2 ) 


90% 

99% 


Table 5-1 | 

Point 

Stress 

Temperature 

Joint Probability of Exceeding 
Stress and Temperature 

1 

p f 

T i 

Rj = 0.01 

2 

p D 

T i 

r 2 = 0.001 

3 

P D 

T 

u 

r 3 = 0.0001 

4 

P i 

T 

u 

r 4 = 0.001 


From the example of Table 5-1 and the joint probability equations it can 
be seen that the joint probabilities of exceeding points 2 and 3 are less than the 
probability of exceeding the design stress. 

For the limit conditions of stress and temperature, point 1, 


R ! - R(p D ) 


If 


R( Pi ) • Rdy >R(p d ) 


In summary, the above discussion shows that if 


R (Tjg) ^ R( Pj ?) 

and 

R(T u ) * R(p D ) 
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the joint probabilities of exceeding either point 2 or point 4 is less than the 
probability of exceeding the ultimate load. Also, the joint probability of 
exceeding both the ultimate load and the ultimate temperature (point 3) is 
much smaller than the probability of exceeding the ultimate load, and design 
for the conditions of point 4 would be excessively conservative. 

5.3.2 Design for Ultimate Load Failure 

The reliability discussion of the previous section suggests an approach 
to design for ultimate load failure at elevated temperature. This approach 
considers failure of the damaged and the undamaged structure and also con- 
siders the sensitivity of the structural integrity to temperature and service 
life. 


In discussing the approach it is assumed that the following criteria are 
known. 


1. The limit loads from which component stresses and loads 
can be computed. 

2. The thermal environment from which component tempera- 
tures can be computed. 

3. Design factors of safety for the mechanical loads from 
which: 


D 


= K 


D **1 


Kq can be either yield or ultimate factors, for example: 


K =1.40 
u 

K = 1.10 

y 

4. The expected service life, t , or the number of flights or 
load cycles, N q , where ° 

t = N /w 

o o' 

and u is the loading frequency in flights per unit time. 
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5. A design life factor, k^, which gives a design life of 



t 

o 


6. A computed temperature, T^, on the component which has 

a probability of occurrence approximately the same as that 
of the limit load, p^. 

7. An ultimate temperature, T , which has a very low prob- 
ability of occurrence. 


T = T . + AT 
u i 


where AT is established after considering trajectory 
considerations, the heating environment, and accuracy 
of computing the resulting component temperature 
spectrum. 


The proposed approach is as follows: 


Step 1 — Ultimate Load Failure of the Undamaged Structure 


Figure 5-lla represents a plane in the failure surface of Fig. 5-10a 
corresponding to the first flight, 



= 0 


a. For the limit temperature, T^, compute the allowable 
load, point 5, 


p A / P| (T i ) * 


Compare the allowable load and the design load (point 2), 
and compute the load margin of safety, M.S. 


M.S. = P 5 /P 2 - 1 


The M.S. should be positive. 
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b. For the limit load compute the allowable temperature, 

T A (P A /Pi=l> = T 6 

Compare the allowable and ultimate temperatures and 
compute the temperature margin (T.M.) 

T.M. = T 6 - T 4 

The T.M. should be positive. 

Alternatively, if it is easier to compute the allowable load at T , a load 
M.S. can be computed for conditions of ultimate temperature and limit load: 

M.S. = P 7 /P 4 - 1 

If all margins are positive the undamaged structure should survive the 
ultimate conditions. Observe that the approach does not require that the com- 
plete failure surface be known. Only two points (5 and 6 or 5 and 7) are required. 

Step 2 — Ultimate Failure at t/t Q = 1 


Figure 5- lib is similar to t-lla except that the failure surface repre- 
sents the allowable load for a service life factor, 

t/t = 1 

and Kj-j and have been replaced by and T^, respectively. is a F.S. 

which defines the design load and Tj is the design temperature for the service 
life factor of one. 

K 1< K D 

and 

T < T. < T 
Si— 1 — u 
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or 

Tj = AT 


The factors Kj and kj^ (or T^) should be specified in the design criteria. 
A preliminary recommendation is to let 


when 

and 


K 1 

= 1.2 

K 

= 1.4 

u 

k l 

= 0.5 


The load M.S. and T.M. are computed in the same manner as for the 
undamaged structure. That is 


M.S. = P 5 /P 2 - 1 

and 

T.M. = T 6 - T 4 

Some hidden conservatism is included in the above margins. The calcu- 
lations for the failure surface, after accumulating damage corresponding to a 
service life factor of one, are based on the factored limit loads and temperature, 
and the assumption is made that these loads are attained on each flight. In 
reality limit loads are the maximum expected loads in normal service and most 
flights will be subjected to loads less than the limit loads. Therefore, the ac- 
cumulated damage should be less than predicted. 


Step 3 — Determine Service Life Margin 

Figure 5- 11c represents the failure surface at 



the design life factor. 
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The scatter in the service life of structures subjected to cyclic and sus- 
tained loads at elevated temperature is notoriously large. It is desired to 
determine the allowable life for limit load and temperature conditions: 

t A^ t o T i ) 

If the allowable life is determined, a service life margin (S.L.M.) can be 
defined 

S.L.M. = t j i / t 0 * k D 

The S.L.M. should be positive. 

If it is more convenient the existance of a positive S.L.M. can be verified 
by computing a positive load M.S. for limit load and temperature. On Fig. 5- 11c 
this is shown by the relationship of points 1 and 5: 

M.S. = p 5 /pj - 1 

5.4 DESIGN FOR PERMANENT DEFORMATIONS 

Permanent deformations from creep or plastic strain can result in 
functional failure prior to rupture or ultimate load failure. Methods of 
predicting creep resulting from sustained loads at elevated temperature are 
discussed in Section 2 and a creep analysis program is described in Appendix B. 

The difficulties of accurately predicting creep deformations increase 
when a structure must be designed for very long life and very low creep rates. 
Most creep data are for relatively high creep rates and the test times are too 
short to be relevant. Utilization of existing data requires extrapolating the 
data for both creep rate and exposure times. As the required extrapolation 
can approach two orders of magnitude, confidence in the extrapolated values 
becomes small. As an example, the design goal for the American SST was 
50,000 flight hours (Ref. 5-9), whereas most creep tests are limited to about 
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500 hours. Additional uncertainties are introduced when the primary creep 
range is considered as little data are available and the data exhibit considerable 
scatter. 

Freudenthal (Ref. 5-2) considers the primary creep range of principal 
importance in contributing to the permissible creep but generally unimportant 
in contributing to damage. The relative importance of primary and secondary 
creep on allowable deformations also vary with the stress and temperature 
conditions. 

5.4.1 Creep Deformation Surfaces 

The dependence of creep deformations on the cyclic stresses is generally 
small (Ref. 5-3 ) as illustrated for M252 alloy at 1350F (Figs. 5-3 and 5-4). 

The design allowable creep condition can be considered independently of the 
ultimate load failure criteria of Section 5.3. 

Figures 5-12 through 5-18 represent typical creep data from MIL-HDBK- 
5B. Figure 5-12 represents creep data in the form of a nomogram for Rene 41 
alloy sheet. Figures 5-13 and 5-14 show short term, stress-time creep data 
for 7075 aluminum, and Figs. 5-15 through 5-18 are long term, stress-time 
creep data for 7075 aluminum. 

Figure 5-19 shows constant creep surfaces for 7075 aluminum which 
were constructed from the data of Figs. 5-13 through 5-18. The constant 
creep surfaces are shown on a stress-temperature-time plot. 

Note that, for high temperatures and long exposure times, the creep 
strains are very sensitive to small variations in either the temperature or 
stress environment. This observation is consistant with the analysis of 
Section 2.4, which considered the case of a small random temperature fluc- 
tuation about the mean temperature, and example problem 5 in Section 2.8. 

The example problem showed that the predicted time for a given creep in a 
titanium beam was reduced from 1750 to 673 hours when a random temperature 
fluctuation of +25°F was added to a 750°F mean temperature. 
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FIGURE 6.3.8. 1 .7. Typical creep properties of Rent 41 alloy sheet. 


Fig. 5-12 - Typical MIL-HDBK-5B Creep Data 
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FIGURE 3.7.2.1.7(b). Creep data for 7075-T6 aluminum alloy (clad sheet) at 400 F. 
Deformation includes thermal expansion of 0.43 percent. Heating rate 75 F per second. 

Fig. 5-13 - Typical MIL-HDBK-5B Creep Data 
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FIGURE 3.7.2.1.7(c). Creep data for 7075-T6 aluminum alloy (clad sheet) at 500 F. 
Deformation includes thermal expansion of 0. 63 percent. Heating rate 75 to 100 F per second 
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FIGURE 3.7.2.1.7(d). Creep data for 7075-T6 aluminum alloy (clad sheet) at 600 F. 
Deformation includes thermal expansion of 0. 74 percent. Heating rate 80 to 90 F per second. 

Fig. 5-14 - Typical MIL-HDBK-5B Creep Data 
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FIGURE 3.7.2.1.7(e). Creep and stress-rupture properties of wrought 707S-T6 aluminum alloy at 94 F. 


Fig. 5-15 - Typical MIL-HDBK-5B Creep Data 
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FIGURE 3.7.2.1.7(f)- Creep and stress-rupture properties of wrought 7075-T6 aluminum alloy at 211 F. 
Fig. 5-16 - Typical MIL-HDBK-5B Creep Data 
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FIGURE 3.7.2.1.7(g). Creep and stress-rupture properties of wrought 7075-T6 aluminum alloy at 300 F. 
Fig. 5-17 - Typical MIL-HDBK-5B Creep Data 
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FIGURE 3.7.2.1.7(h). Creep and stress-rupture properties of wrought 7075-T76 aluminum alloy at 375 F. 


Fig. 5-18 - Typical MIL-HDBK-5B Creep Data 
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5.4.2 Design Factors for Allowable Creep 

The proposed approach to design for limiting creep deformations is 
similar to the approach to design for ultimate load failure. The allowable 
creep design factors are used to assure that adequate margins on load, tem- 
perature and time exist between the design conditions and the allowable creep 
surface. 

Step 1 — Creep at a Service Life Factor of One 

At a time equal to a service life factor of one, determine the allowable 
creep stress (the stress which results in the allowable creep strain) in the 
component at the limit temperature. 


F A (T r l / l o = » 


Compute the design stress, f^, 


D 


= K 


where 


K c is a creep F.S. 

f^ is a stress corresponding to limit load. 


Compute the creep stress M.S. 


F A (T i' l /*o = ** , 

M.S. = 7 - 1 

D 


Next, determine the allowable temperature which corresponds to f 


f 


t a VAo =1 > 
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Compute a creep temperature margin 

T.M. = T a (f r t/t o = 1) - T £ 

The T.M. should be greater than a temperature safety margin, AT c> which is 
analogous to a F.S. on mechanical loads. 


Step 2 — Service Lifetime to Limiting Creep 


The final check consists of computing a life margin(S.L.M.)cor responding 
to conditions of limit temperature and stress. Compute the allowable service 
life factor. 


‘aA„ <VV 


Compute the S.L.M. from 


S.L.M. - t A /t o (f|, T^) - t D A Q 

where t^/t Q is the design service life factor. The S.L.M. should be positive. 
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Appendix A 

CREEP DATA FOR METALLIC MATERIALS 


ft-f (*■) 
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Appendix A 

Material data for the group of metals listed in Table A-l have been 
gathered from references listed in Section A.l for each of the materials. 

One of the most striking features of these graphs, shown in Section A. 2 is 
the extreme senstivity of some of the materials to variation of the data with 
stress and temperature. 

A.l MATERIAL REFERENCES 

At elevated temperatures, material properties become greatly tempera- 
ture and time dependent. These temperature-time-load characteristics must 
be known with a fair degree of confidence in the accuracy and reliability for 
the data to be useful in design work. In design handbooks such as MIL-HDBK-5B 
the following material properties can be found for most aerospace materials: 

• Thermal conductivity, thermal expansion coefficient and specific heat 
versus temperature 

• Strength versus temperature for various exposure times 

• Charts for strength at test temperatures versus exposure time and 
exposure temperature 

• Strength versus time for given total deformation (creep) at constant 
temperature 

• Constant life diagrams for various exposure times at constant 
temperature 

• Stress -strain diagrams for various temperatures 

• Elastic modulus versus temperature 

• Tangent moduli versus stress for given temperatures 

• Nomogram for computing creep strain for given exposure time to 
given stress and temperature. 
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High temperature creep is one area where sufficient data are lacking 
for quite a few structural materials and in particular those materials applicable 
to high temperature ranges. Table A-l. 

During the course of this study several creep references have been com- 
piled and have been listed here as additional sources of information to supple- 
ment MIL-HDBK-5B. The materials are listed alphabetically by their common 
name with a short description of the information in each reference. 

Aluminum 


2024-T3 

Heimerl, George J., and Arthur J. McEvily, Jr., "Generalized Master Curves 
for Creep and Rupture," NACA TN 4112,, October 1957. 

Master curves for creep rupture and minimum creep rate. 


2219 


Aerospace Structural Metals Handbook, Air Force Materials Laboratory, 
Wright Patterson AFB, Ohio, 1970. 

cr versus t for percent creep and rupture for various temperatures. 

Alloy Digest , Engineering Alloys Digest, Inc., Upper Montclair, N. J. 

Data table cr versus t, percent creep in specific time. 

7075-T6 

Heimerl, George J., and Arthur J. McEvily, Jr., "Generalized Master Curves 
for Creep and Rupture," NACA TN 4112, October 1957. 

Master curves for creep rupture and minimum creep rate. 


C olumbium (Niobium) 


CB-753 


Conway, J. B., and P. N. Flagella, "Creep -Rupture Data for the Refractory 
Metals to High Temperatures," Gordon and Breach Science Publishers, 1971. 

a versus t for various percent creep for two temperatures 
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Table A -1 

TEMPERATURE AND MATERIALS 
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Cb-752 

Conway, J. B., and P. N. Flagella, "Creep -Rupture Data for the Refractory 
Metals to High Temperatures, " Gordon and Breach Science Publishers, 1971. 

a versus t for percent creep for three temperatures 


D-43 

Conway, J. B., and P. N. Flagella, "Creep-Rupture Data for the Refractory 
Metals to High Temperatures, " Gordon and Breach Science Publishers, 1971. 

a versus t for percent creep at three temperatures 

Aerospace Structural Metals Handbook, Air Force Materials Laboratory, 
Wright-Patterson AFB, Ohio, 1970. 

Larson-Miller parameter and creep rupture curves for three temperatures. 


Inconel (Nickel-base alloy) 


Inconel 702 


Aerospace Structural Metals Handbook, Air Force Materials Laboratory, 
Wright-Patterson AFB, Ohio, 1970. 

Larson-Miller Parameter of 0.2% creep and creep rupture, 1970. 

''Research Investigation to Determine Mechanical Properties of Nickel and Cobalt- 
Base Alloys for Inclusion in Military Handbook-5,* ML-TDR-64-116, Air Force 
Materials Laboratory, Wright-Patterson Air Force Base, Ohio, October 1964. 

a versus t at various T for various percent deformations 

Inconel X 

Heimerl, George J., and Arthur J. McEvily, Jr., "Generalized Master Curves 
for Creep and Rupture," NACA TN 4112, October 1957. 

Master curves for creep rupture and minimum creep rate. 

L605 (Cobalt base alloy) 

"Joint International Conference on Creep, 1963," sponsored by ASME, ASTM 
and the Institution of Mechanical Engineers, 1963. 

a versus t for 0.5% creep for three T 
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Aerospace Structural Metals Handbook, Air Force Materials Laboratory, 
Wright-Patterson AFB, Ohio, 1970. 

Manson-Haferd parameter for 0.5% and 1.0% plastic strain 

a versus percent plastic e at 1800°F for various t intervals 

Alloy Digest, Engineering Alloys Digest, Inc., Upper Montclair, N. J. 

^Research Investigation to Determine Mechanical Properties of Nickel and 
Cobalt-Base Alloys for Inclusion in Military Handbook-5? ML-TDR-64-116, 
Air Force Materials Laboratory, Wright-Patterson Air Force Base, Ohio, 
October 1964. 

a versus t at various T for various percent deformation 


Molybdenum 


TZC Alloy 

Sawyer, J. C., and E. A. Steigerwald, "Generation of Long Time Creep Data 
of Refractory Alloys at Elevated Temperatures," TRW, Inc., Cleveland, Ohio, 
(NASA CR-1115). 

Larson-Miller parameter, 0.5% creep 


TZM Alloy 

Aerospace Structural Metals Handbook, Air Force Materials Laboratory, 
Wright-Patterson AFB, Ohio, 1970. 

Creep rupture curves a versus t for Ts. 

o versus t for percent creep at 2000°F 


Rene 1 41 (Nickel-base alloy) 

"Joint International Conference on Creep, 1963," sponsored by ASME, ASTM 
and the Institution of Mechanical Engineers, 1963. 

a versus e for various T 

Aerospace Structural Materials Handbook, Air Force Materials Laboratory, 
Wright-Patterson AFB, Ohio, 1970. 

Manson-Succop Parameter - Creep Rupture, 1970 edition. 

Larson-Miller Parameter - 2% creep + a versus t for various T, 

1966 data. 
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Gluck, J.V., and James W. Freeman, "Effect of Creep-Exposure on Mechani- 
cal Properties of Rene' 41, " ASD TR 61-73, Air Force Material Laboratory, 
Wright-Patterson Air Force Base, Ohio, August 1961. 

Plots of e versus t - obtain strain intercept data 

"Research Investigation to Determine Mechanical Properties of Nickel and 
Cobalt-Base Alloys for Inclusion in Military Handbook-5, " ML-TDR-64- 1 16, 
Air Force Materials Laboratory, Wright-Patterson Air Force Base, Ohio, 
October 1964. 

a versus t at various T for various percent deformation 


Steel 


Carbon-Molybdenum Steel and 18-8 Cb Stainless Steel 

Heimerl, George J., and Arthur J. McEvily, Jr., "Generalized Master Curves 
for Creep and Rupture," NACA TN 4112, October 1957. 

Master curves for creep rupture and minimum creep rate. 


T antalum 


T-lll Alloy, and T-222 

Sawyer, J. C., and E. A. Steigerwald, "Generation of Long Time Creep Data of 
Refractory Alloys at Elevated Temperatures," TRW, Inc., Cleveland, Ohio, 
(NASA CR-1115). 

Larson-Miller and Marson-Haferd Parameters of 1% creep at 3000°F 

Plots of percent creep versus t for different T and a (T-lll) 

Aerospace Structural Metals Handbook, Air Force Materials Laboratory, 
Wright-Patterson AFB, Ohio, 1970. 

Stress rupture curves at 2400°F 

TD-Ni Cr 


Hirschberg, M. H., David A. Spera and Stanley J. Klima, "Cyclic Creep and 
Fatigue of TD-Ni Cr, TD-Ni, and Ni Cr Sheet at 1200C, " NASA TD D-6649, 
February 1972. 

Creep rupture curves (1200°C) 
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Titanium 


Ti-6Al-4V 

"Joint International Conference on Creep, 1963, " sponsored by ASME, ASTM 
and the Institution of Mechanical Engineers, 1963. 

e versus t for various T at constant a 

Alloy Digest, Engineering Alloys Digest, Inc., Upper Montclair, N. J. 

Plot of total creep versus t for two T and two a. 

Aerospace Structural Metals Handbook, Air Force Materials Laboratory, 
Wright -Patters on AFB, Ohio, 1970. 

a versus t for various T at specific percent creep. 

Tungsten 


Conway, J. B., and P. N. Flagella, "Creep-Rupture Data for the Refractory 
Metals to High Temperatures, " Gordon and Breach Science Publishers, 1971. 

e versus t curves 

Creep rupture curves 

a versus e for various temperatures 

Larson-Miller parameter for rupture 

Aerospace Structural Metals Handbook, Air Force Materials Laboratory, 
Wright -Patters on AFB, Ohio, 1970. 

Creep rupture curves for four temperatures. 


A. 2 CREEP CURVES FOR CONSTANT STRESS AND TEMPERATURE 


In this section curves of the type of Fig. 2-2 are presented for a variety 
of metallic materials. This is the form in which the data can be used for creep 
calculations (Figs. A- 1 through A -9). 
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Fig.A-5 - Inconel 702 
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Fig.A-7 - Molybdenum TZC Alloy 


LMSC-HREC TR D306579 








LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D306579 






























LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D306579 

















Appendix B 

CREEP ANALYSIS PROGRAM 


13 ' * O) 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D306579 


Appendix B 

The program input data, sample problem, program listing and the output 
for the sample problem are included in this appendix. 

B. 1 PROGRAM INPUT 

The required data input and format to be followed are outlined in the fol- 
lowing statements. Note, at least two cards are required in several sections 
of data. Cards 1,2,5 and 6. Number of cards depend on number of data sets 
to be entered. 

Cards 1 


NSIG, DTITLE, (SIGMA(I), P(I), 1=1, NSIG) 
FORMAT (15, 10A6/(8E10.4)) 


NSIG 
DTITLE 
SIGMA (I) 

P(I) 


Number of sets of data points to be entered 
Title or label for data 

Stress value from master creep rupture parameter 
curve (psi), largest stress value first. 

Corresponding parameter value 


*NOTE: Correct parameter equation must be input into the program correspond- 
ing to the parameter data input, see Section 2.9. 


Cards 2 


NC, STRN, CTITLE, (SIGCR(I), PCR(I), 1=1, NC) 
FORMAT (15, E10.4, 9A6/(8E10.4) 

NC Number of sets of data points to be entered 
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STRN 

CTITLE 

SIGCR(I) 

PCR(I) 


Percent creep strain for which the parameter data 
applies, i.e., 0.2% creep strain, etc. 

Title or label for data 

Stress value from master creep parameter curve 
(psi), largest stress value first 

Corresponding parameter value 


*See note in Cards 1 


Card 3 

STRMAX 
FORMAT (E10.4) 


STRMAX Maximum strain value allowed (in. /in. ), card can 

left blank. 


Card 4 

LIREQ 

FORMAT (E10.4) 

LIREQ Total life required (hrs) 


Cards 5 


NS, TITLE, (STRESS(J), TEMP(J), TIME(J), J=l, NS) 
FORMAT (15, 10A6/(3E10.4)) 


NS 

TITLE 
■ STRESS(J) 
TEMP(J) 
TIME(J) 


Number of load cases 
Descriptive label for load cases 
Stress level (psi) 

Temperature (°F) 

Time period stress and temperature applied (hrs) 


*NOTE: A stress, temperature and time constitute a load case. 
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Card s 6 


NSTR, VSTRES(I), VTEMP(I), (VTIM(J, I), VSR(J, I), J=l, NSTR) 
FORMAT (15, 5X, E10.4, E10.4/(8E10.4)) 


NSTR 
VSTRES(I) 
VTEMP(I) 
VTIM(J, I) 
VSR(J, I) 


Number of data sets 

Stress level for which data applies (psi) 

Temperature level for which data applies (°F) 

Time values from plot of transient strain vs time (hr) 
Corresponding transient strain value (in. /in. ) 


*NOTE: NS sets of transient creep data must be input, 1=1, NS. The order of 
the NS sets of data does not have to be the same as the order in the 
load cases. The program will automatically search for the correct 
data based on VSTRES and VTEMP. 


B.2 SAMPLE PROBLEM 


An example problem has been constructed using Rene' 41 material and a 
variable loading condition. The following stress history was input. 

Table 1 

LOAD CONDITIONS 


NS 

Stress 

Temperature 

Time 


(psi) 

(°F) 

(br) 

1 

25,000 

1600 

4.0 

2 

35,000 

1500 

3.5 

3 

85,000 

1300 

4.5 

4 

100,000 

1300 

1.5 


The Manson-Succop creep rupture parameter for Rene' 41 was taken from 
Ref. B-l and the Larson-Miller creep parameter for 0.2% creep strain from 
Ref. B-2. A maximum allowable strain of 0.15 in. /in. and a required life of 
200 hours were arbitrarily chosen. 
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Transient creep strain data were determined for Rene’ 41 from the data 
presented in Ref. B-3. These material data, Fig. B-l, were derived only for 
use in writing and checking out the operation of the digital computer program 
and should be used with caution. 

The output data and plots for this problem are given in Section B.4. 
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SECTION B.3 PROGRAM LISTING 


DATE 010373 PAGE ? 


OOUl 


002270 

920L 

0000 

I 

007055 

bcox 

0000 

R 

007071 

6CDY 

0000 

R 

007413 

COEFFM 

oooo 

R 

0U1666 

CRATE 

00U0 

K 

007373 

CRATE1 

0000 

R 

007401 

CRATL2 

0000 

R 

001750.., CREEP 

0000 

R 

007357 

CREEPT 

0000 

R 

002222 

CTITLE 

OOdO 

H 

007605 

CX 

OOOO 

R 

007406 

DELSTR 

0000 

R 

0022.11) 

dtitlf 

0000 

R 

007131 

EXPN 

oooo 

R 

007414 

FS 

OOuO 

K 

000764 

FSUM 

0000 

R 

000702 

PT 

0000 

R 

007«ift4 

FX 

0000 

I 

007347 

I 

oooo 

I 

007354 

J 

OOuO 

I 

007361 

K 

0000 

0 

007155 

M 

0000 

O' 

007225 

K2 

0000 

I 

007356 

L 

oooo 

R 

007345 

LIREQ 

OOUO 

1 

007362 

N 

0000 

I 

00735U 

NC 

0000 

I 

007415 

NP 

0000 

I 

007353 

NS 

oooo 

1 

0U7346 

NSIG 

OOuO 

2 

007355 

NSTR 

uooo 

I 

003301 

NV 

00 ,Djft 

K Qooofva 

P 

0000 

R 

001440 

PC 

oooo 

R 

001274 

PCR 

OOUO 

H. 

007371 

Prtli 

0000 

R 

007377 

PRTL2 

onoo 

R 

007372 

PPTl 

0000 

R 

007400 

PMT2 

oooo 

K 

000454 

PX 

OOUO 

H 

007370 

PI 

0000 

R 

007376 

P2 

0000 

D 

007275 

G 

oouo 

R 

001356 

sigclg 

oooo 

R 

001212 

SIGCR 

OOuO tR 

000144 

SIGLOG 

occo 

R 

000000 

SIGMA 

0000 

R 

001130 

stlog 

0000 

R 

007375 

STLl 

oooo 

R 

002233 

STPRI 

OOUO 

H 

007363 

S7PHI1 

0000 

R 

007364 

STPRI2 

0000 

R 

007402 

str 

0000 

R 

000226 

STRESS 

oooo 

R 

007352 

STRMAX 

OOuO 

H 

007351 

5TRN 

0000 

R 

00736b 

STRO 

0000 

R 

007374 

ST1 

0000 

R 

006627 

SY 

oooo 

R 

OU 1 604 

TC 

OOUO 

K 

001522 

1 CLOG 

0000 

R 

000310 

TEMP 

0000 

R 

000372 

time 

0000 

R 

005725 

TIME1 

onoo 

R 

006153 

TIME2 

OOuO 

K 

007410 

TISUM 

uooo 

R 

002176 

TITLE 

0000 

R 

000536 

TLOG 

0000 

R 

007412 

TMAX 

oooo 

R 

007407 

TO 

OOuO 

R 

000620 

TR 

0000 

R 

007105 

TREM 

0000 

R 

007360 

TS 

0000 

R 

002461 

TSTR 

oooo 

R 

006401 

TX 

OOuO 

R 

00 74 03 

1XX 

0000 

R 

002707 

12 

0000 

R 

004544 

VSR 

0000 

R 

003135 

VSTRES 

oooo 

R 

003217 

VTEMP 

OOUO 

R 

003363 

VTIM 

0000 

R 

007411 

XLIFE 

uooo 

R 

001046 

XMAR 

0000 

R 

007366 

XTl 

oooo 

R 

007367 

XT2 


ooxui 

I** 



DIMENSION SIGMA (50) «P(SO) 

0J1U3 

2* 



DIMENSION SIGL0G(50)»STRESS(50> »TENP(50>#TIME<50)#PX(50>*TLOG<S0) 

OU104 

3“ 



DIMENSION TR (50) #PT (50) ,FSUM(50> »XMAR ("50) » STLOG (50) 

0U105 

4 * 



DIMENSION SIGCk ( 50 ) »PCH ( 50 ) > SIGCLG (50 ) » PC<50>» TCLOG (50 ) , TC (50) » 

UUJ05 
'f 00106 

5* 



ICR ATE ( 50) » CREEP ( 5U» 3) * I ITLE ( 10 ) *DTITLE (10) , CTITLE <9 > 

6* 



DIMENSION 5TPRI(5U»3> » TSTR ( 50 # 3 > * T2 (50 # 3) 

o 0010? 

7* 



DIMENSION VSTRES (50) » VTEMP (50) ,NV(50) » VTIM(25.25) » VSR (25*25) * 

00107 




1TIME1(50»3)»TIME2(50*3) 

00110 

9* 



DIMENSION TX (50,3) >SY (50,3) *BCDX ( 12 ) » BCDY ( 12 ) 

UOlll 

10* 



DIMENSION TRtM(20) * EXPN (20) #K1 ( 20) »K2 (20 ) *0 (20) 

00112 

11* 



REAL LIRtQ 

OU 1 13 

12* 



INTEGER bCUX 

0,U14 , 

13 * 



DOUBLE PRECISION K1»K2,Q 

UU114 

14* 

C 


THE FOLLOWING DATA ARE THE MANSON-SUCCOP PARAMETER AND STRESS 

0 M 14 

15* 

c 


VALUtS FOR RENE Al HAR-CREEP RUPTURE (AERO. STRUCT. METALS HDBK-1970) 

00114 

16* 

c 


P = LOG ( TIME ) + 0.0106 TEMP (DEG, F) 

0U114 

17* 

c 


THE FIRST STRESS VALUE IS THE LARGEST STRESS VALUE IN THE TABLE 

00114 

)8* 

c 



OU 1 14 

19* 

c 


THE CPEEP RUPTURE DATA ARE READ IN SIGMA AND P. SIGMA= STRESS AND 

Oil 114 

20* 

c 


p= PARAMETER VALUE. FROM MASTER CURVE 

OUJ 14 

21* 

c 



00115 

22** 



READ(5»lU0> NSIGrUlITLE' (SIGMA ( I ) ,P ( I ) * 1 = 1 , NSIG) 

00131 

23* 


ioo 

FORMAT ( 15 » 10A6/ (6E.10.A ) > 

U 0132 

24* 



IF (NSIO.EQ . 0 ) GO 10 205 

OU 134 

25* 



WRITE (6. 101 ) DIITLt 

00142 

26* 


101 

F0nMAT(1h1»10Ao) 

0U143 

27“ 



WRITE (6. 105) ( (SIGMA (I) , P( I )) , 1=1 , NSIG) 

00162 

2H* 


105 

FORMAT (liiO,»El5.6> 

0U153 

29* 



co no r=i.f;siG 

0 0166 

30- 



SlGLOGl I)=AL0G10(S16MA(I) > 

00157 

31 ♦ 


110 

CONTINUE 

0 ul61 

32’*' 


205 

CONT INUE 

0 u 161 

33* 

L 


. 


34 * 

C 


THE CREEP MASTtR CURVE DATA ARE READ IN SIGCR AND PCR, STKN= PER- 

o j i b 1 

3b* 

C 


CI.ImI CREEP POP. WHICH THE DATA APPLIES , SIGCR= STRESS And PCR= PA- 

00161 

36 * 

u 


RA., Elf H vALUt. FROM THE MASTER CURVE 
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00161 

37* 

C 



0U162 

38* 



REAO (5 > 210 ) NC.STKN.CTITLE* (SIGCR(I).PCR(I) >1=1, NC) 

00177 

39* 


210 

FORMAT ( IS .E10.4.9A6/(8E10.4>) 

00200 

40* 



IF(NC.EQ.O) GO TO 240 

00202 

41* 



WRITE(6»220> CTITLE.STRN 

00211 

42* 


220 

FORMAT ( 1H0 . 9A6/ 1H0 » 23HAM0UNT OF CREEP STRAINS. E15.8> 

00212 

43* 



WRITE (6.230 ) (SIGLR(I) »PCR(I> .J^rl^Nq,), i 

00221 

44* 


230 

FORMAT (1H0. BE lb . 6 ) 

0 U 222 

>5* 



00 2S0 1=1. NC 

00225 

, 4b* 



SIbCLG( I )=ALOG10(bIGCR( I ) ) 

00226 

47* 


250 

CONTINUE 

0 0230 

48* 


240 

CONTINUE 

0U231 

49* 



REAU(5r2o0) STRMAa 

00234 

50* 


260 

FORMAT (E10 . 4 ) 

00235 

♦DIAGNOSTIC* 

THE TEST FOR EQUALITY BETWEEN NON-INTEGERS MAY NOT BE MEANINGFUL. 

00235 

bl* 



IF(STRMAX.EO.O.O) GO TO 270 

0U237 

b 2 * 



WRITE (b.*65> STRMAX 

00242 

53* 


265 

FORMAT (1H0.22HMAX STRAIN ALLOWABLES. E15. 8) 

UQ243 

54* 


270 

CONTINUE 

0U244 

b5* 



READ l 5 .405) LIREQ 

0*i247 

56* 


405 

FORMAT ( E10.4) 

00250 

b7" 



WRITE (6»40b) LIREQ 

00253 

58* 


406 

FORMAT (lH 0 . 14 HLIFt REQUlRE0=.E15.a) 

00253 

69* 

c 



00253 

oO* 

c 


THE FOLLOWING VALUES ARE THE APPLIED LOADS 

002*53 

bl * 

c 


STRESS IS THE APPLIED STRESS LEVEL AT TEMP. AND THE PERlOO OF TIME 

fcj 00253 

b 2 * 

c 


THAT THE STRESS IS APPLIED . STRESSIPSI ) .TEMP(F) . TIME(HR) 

* 00253 

63* 

c 



0U254 

64* 



READ (5» 115) NS. TITLE. (STRESS CJ) .TEMP(J) .TIME(J) .J=1»NS) 

002/1 

65* 


115 

FORMAT (I5*10A6/(3LlO. 4) ) 

0 U 2/2 

06 * 



WR1TL ( 6 . 120 ) 

002/4 

67* 


120 

FORMAT (1H1. 54h THt FOLLOWING VALUES ARE THE PRESCRIBED LOAOING CY 

0 02 74 

08 * 



lCLE/JX.1 1HSTRESS (RSI ) . 11X .7HTEMP (F ) . 14X » BHTIME (HR) ) 

0U275 

j o9* 



WR1TL ( 6 » 125 ) (STRtSS(J) . TEMP ( J) . T IME ( J) .J=1»NS) 

00 505 

70* 


125 

FORMAT (1H0. 3 (E15.B.5X)) 

0030b 

71* 



WRITL (6.501 ) 

00310 

72* 


501 

FORMAT (1H1. 99HTHfc FOLLOWING DATA ARE THE TRANSIENT CREEP STRAIN V 

0 U 5 1 0 

73* 


2ALUES FOR THE LOAD STRESS AND TEMPERATURE VALUES) 

00311 

74* 



DO 500 1=1. NS 

00314 

75* 



READ(5»5U5)NSTR,VSTRES(I).VTEMP(1)» ( VTIM ( J. I ) . VSR ( J » I ) . J=1 .NSTR ) 

00326 

76* 


505 

FORMAT ( 15. 5X.E10.4 .E10 .4/ (8E10.4 ) 1 

0 0 327 

n * 



WRIT) (b.50b)NS!R» v5TRFS« 1 ) » VTFMP( I ) . (VTTM(J.I) » VSR( J.I ) .0*1 .NSTR) 

00341 

ih* 


506 

FORMAT (1H0.3HNV=.15»5X»7HS]RESS=»E15.8»5X»SH1EMP=,E15.B/(6(E15.8* 

00541 

79* 



1 5X) ) ) 

00342 

60* 



NV(I)=NSTR 

0U343 

bl* 


500 

CONTINUE 

00345 

o2* 



DO 131 L=l»3 

0U3SU 

b3* 



CREEPT=0 . 0 

0U351 

d4* 



TS=0.0 

00352 

6b* 



00 130 1=1. NS 

00355 

bb* 



K = 0 

O'Obb 

b7 * 



IF (NSlb.LG.O) GO TO 70 

00360 

68* 



IF (SI RLS5 ( I ) -SIGMA ( 1 ) ) SO. 50 .60 

il.; 



50 

cnnTlNUF 

»i 4 

VO* 



11 (MRI SS(1)-*..KMA(N‘ J 10) >60.70.70 

ij . 1 5 7 

VI* 


60 

WRIYL(fa.b5) I .S)KtSS(l) 

M 0 5 ‘ .t 

S2 * 


bb 

FORMAT ( 4X . 13HSTRF.5S VALUL(«I2. 2H> = . E 15.8 . 2X » 3PHEXCEEDS THE R*NGE 
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00573 

*3* 


10F THE RUPTURE TABLE) 

0J374 

y 4* 


TR(I)=0.0 

00575 

9b* 


FT ( I )=0.0 

00376 

9 6+ 


K=1 

00377 

97* 

70 

CONTINUE 

00400 

98* 


IF (NC.EQ.O) GO TO 370 

00402 

99* 


IF(STRESS(I)-SIGCR<1) )350#350»,36'0' 

00405 

luO* 

350 

CONTINUE 

00406 

1 Ul + 


IF(STRESS(I)-SIGCK(NC) >360*370*370 

0041 X 

• lu2 + 

360 

WRITE (6*36b) I*STKtSSU> 

0 04 1 5 

103* 

365 

FORMAT <4X»13HSTRESS VALUE< f I2*2H)=*E15.8*2X*36HEXCEEDS THE RANGE 

00415 

104* 


IF THE CREEP TABLE) 

00416 

lu5* 


CREEP ( 1 >L) =0 .0 

00417 

106* 


GO TO 130 

00420 

107* 

370 

CONTINUE 

UU421 

108* 


STLOG<I>=ALOG10<STRE5S(I) > 

00422 

109* 


IF (NSIti.EQ.O) GO 10 310 

0 04 24 

110 + 


IF(K.GT.O) GO TO 320 

00426 

111 + 


CALL GIRKSTLOGU) *PX(I> *SlGLOG*P»NSIG> 

0042b 

112* 

C 


00426 

113* 

c 

THi FC’-LOWlNG EQUATION IS THE MANSON-SUCCOP PARAMETER EXPRESSION 

00426 

1 14 * 

c 


00427 

115 + 


TLOG ( I ) SpX ( I ) — ( , 010fl*TEMP ( 1 ) ) 

0U427 

116* 

c **** 

00427 

117* 

c **** 

t0OU43O 

118* 


7R(I)=10.0**TLOGU) 

* 004 31 

119* 


FT(I)=TIME(I)/TR(I) 

00 0 04 32 

120* 


GO TO 320 

00433 

121- 

310 

TR<I)=G.O 

00434 

lc2* 


FT<I)=0*0 

00435 

lc3* 

320 

CONTINUE 

004 36 

124* 


DO 510 N=1#NS 

01)441 

♦UImGNOSTIC* THE TtST FOR EQUALITY BETWEEN NON-INTEGERS MAY NOT BE MEANINGFUL. 

0 0 4 4 i 

1^5* 


IF (S 1RESS ( 1 ) .EG. VSTRES(N) ) GO TO 520 

0 0 4 m 2 

lc6* 


GO TO blO 

00444 

♦DIAGNOSTIC* THE TtST FOR EQUALITY BETWEEN NON-INTEGERS MAY NOT BE MEANINGFUL. 

0 0444 

127* 

520 

IF(Tt>P(l) .EO.VTfertP(N) ) GO TO 530 

00446 

128* 

510 

CONTINUE 

0 0 4 5 0 

129* 

530 

J-NV < N ) 

00451 

150* 


TIKE1 ( I >L)-TS 

00452 

131* 


TlME2«I.L)=TS+TIMt(I) 

00453 

io2* 


IF (L.EQ.2) GO TO B50 

00455 

1 33 * 


IF (L.EQ.3) GO TO 1000 

00457 

1 54 * 

855 

IF(TS.GT.VTIM<0*N) ) GO TO 540 

0 0 4 c 1 

lob* 


CALL INTtRl < 1S.STPHI 1 * VT1M( 1 *N) • VSR ( 1 *N) » J) 

004t>2 

1 36* 


IF (T1NE2( I *L).bT.VlIM(«J*N)) GO TO 513 

0U464 

137* 


CALL INTtRl (TIME2U.L) *STPRI2*VT1M(1»N> . VSR(1*N) .0) 

004t>5 

l3fl* 


STPRI ( I.L)=STPRl2-bTPRll 

004o6 

139* 


IFU.EO.l) GO 10 b08 

00470 

140“ 


GO TO S09 

0)471 

l4l* 

508 

STH0=STPRI1 

004 72 

142 * 


TSTR(l»L)=STPRl(l.l.) 

00473 

143» 

509 

CONTINUE 

:) -4 74 

144-* 


1! (llf>t'«:(I«L).G1.VTIM(J*N)) GO TO 513 1 

.*04 Mi 

14 5* 


'I SI l< < I * L ) - ! M P ( I - 1 * I I+STPRKI.L) 

■jju n 

3. 4 6 ^ 


GO Tc 7b 


o imm i ni > ni.n ~ ' r i - ■ \ ‘W. *r ^ .i t. 
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00502 

ms* 



IF (TSTR (I-l.L).GT.VSR(J'N') GO TO 540 

00504 

Vny* 



CAUL 1NTER1(TSTRU-1*U 'XTl,VSR<l»N>»VTIMUrN>rJ> 

00505 

lbO* 



XTjiSXTl+TIME ( I ) 

00506 

lbi* 



IF (XT2.GT.VTIM(J»N> ) 60 TO 522 

00510 

)b2* 



CALL INTER1 ( XT2 »STPRI2» VTIM ( 1 rN) 'VSR(lfN) »J> 

00511 

lb3* 



TSTR(I*L)=STPR12 

00512 

154* 



GO TC 75 

00513 

155* 


513 

TRtM(I)=TIME2(I»U-VTIM(J»N) 

00514 

156* 



GO TO 523 

00515 

157* 


5 22 

TRLM<I)=xT2-VTlM(v"N) 

00516 

*158* 


523 

STPRI2=VSR( J»N) 

00517 

159* 



IF(L.EG.2> GO TO bl6 

00521 

160* 



STPHI < I .L>=STPKI2-STPRI1 

UU522 

lol* 



GO TO 550 

00523 

lo2* 


540 

STPRI ( 1 rL)-0 .0 

00524 

163* 



GO TO 550 

00525 

lc>4* 


516 

STPRI ( 1 »U=STPRI2-TSTR(I-1*L) 

0 0 526 

lo5* 


550 

CONTINUE 

0052/ 

lb6» 



IF INC. EG. 0) GO TO 300 

00631 

lt>7* 



CALL GIRKSTLOGU) f PC (1). sigclg.pcr.no 

00531 

l&a* 

c 



0 053 1 

lo9* 

c 


THt FOLLOWING EQUATION IS THE LARSON-MILLER PARAMETER 

0 053 1 

l/o* 

c 



00532 

171* 



TCLOG<I>= < (PC(I) *1000. >/(TEMP(l>*460.) »-20. 

00532 

172* 

c 

***♦ 

H0 00532 

173* 

c 

*♦** 

i 0U533 

174* 



TC(I)=10.0**TCL0Gtl) 

^ 00534 

17b* 



CRATE(I)=STRN/TCU) 

00535 

176* 



IF UREM(I) .GT.0.0) GO TO 531 

00537 

177* 



IF(T1PE 1(I'L1. GT.VTlMj'N)) GO TO 511 

0 0541 

178* 


531 

CREEP(I#L)=CRATE<1)*TREM(I> 

.00542 

1 79* 



GO TO 512 

00543 

160* 


511 

CREEP (l»L)=CRATE(l)*TIME(I) 

00544 

lol* 


512 

CONTINUE 

00545 

' 162* 



TSTR ( I 'L ) sTSTR ( I-1'L 1+STPR1 ( 1 »L 1 ♦CREEP ( I»L) 

0 0546 

163“ 



GO TO 75 

00546 

164* 

c 



00546 

165* 

c 

*** 

PAO-MARIN THEOHY 

0 U 5 4 6 

lt>6* 

c 



00647 

167* 

1000 

CONTINUE 

00550 

168* 



CA| L GIR1(STL0G(I)'P1'SI6CLG'PCR'NC) 

0 0661 

164* 



PNl 1 L 1= < CP1*1000.1/(TEMP(1) *460.1 )-20. 

00652 

lyu* 



PMI1=10.U**PMTL1 

00563 

lyl* 



CRATtl=STRN/PMTl 

00554 

ly2* 



STisSTRESSm+bOOO. 

00555 

iy3* 



STL1=AL0610(ST1> 

00556 

ly4* 



CALL GIRl<STLl*P2»bIGCLG'PCRrNC) 

U 0 5 5 7 

iy5* 



PMTL2=< (P2*1000. >/<TEMpU)i-4b0.d-20. 

OU560 

ly6* 



PMT2=10.U**PNTL? 

00561 

iy7* 



CRATE2=SIRN/PMTP 

00562 

lya* 



EXPN (l)-«L<J6(CKATt.2/CRATLl) /ALOG (ST 1 /STRESS (11) 

00563 

lyy* 



Kid) -VSK ( J »N > / t S 1 KFSS ( I 1 **FXPN < I ) ) 

00664 

200* 



K?(I)=U<ATLl/CbTi«tbS(n**EXPN(l > > 

UOboS 

2ol * 



CALL INTER 1 ( 1 IMF ( 1 ) »STl<»VTlM(i»Nl »VSH<1»N>.J> 

0o6f.6 




IF <1,01.1) 0(i TO 

0 0670 

2u3* 



XTim.o 

00571 

2C4* 



XU.=T Il-Edl 
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Q0572 

205* 


GO TO 2999 

00573 

2 ‘ 

2990 

CALL INT£RHTSTRU-l»U*XTl,VSRU*N)»VTIM..<J»N>»J) 

00574 

2u7* 


XT2SXT l+TIMEd ) 

00575 

20ft* 

2999 

CONTINUE 

00576 

209* 


<3(I)=(1./TIME<I) )*AL0G(1./(1.-(STR/(K1(I)*STRESS(I)**EXPN(I))))) 

00577 

210* 


IF (l.EQ.l) GO TO 3150 

OObOl 

211* 


IF (TSTR(I-1»L) ,Gfc.VSR(J*N) > GtQ .'TU «MV»0 

U0o03 

212* 

3100 

IF(XT2.G1 ,VTIM(J»N) ) GO TO 3300 

00605 

213* 


IF(XT2.Lt.VTIM(J»N) ) GO TO 3500 

00607 

<214* 


GO TO 3200 

00610 

215* 

3500 

K2 ( I )-0 • 0 

00b 1 1 

216* 


GO 10 3200 

00612 

217* 

3300 

TXX=VTIM(J»N)-XT1 

00613 

218* 


TREM(I)=XT2-VTIM(u»N> 

0 U o 1 4 

219* 


TSTR(I»L)= (K1(I)*(EXP(-0(I)*XT1)-EXP(-Q(I)*VTIM(J»N)>) +K2(D* 

00o 14 

220* 


l TREM(I) )*( (ABStSlKESS(I) )**(EXPN(I)-1.) )*STRESS<I))+TSTR<I-1.L) 

00615 

221* 


GO TO 3400 

00O16 

222* 

4000 

GU)=0.0 

00617 

223* 


GO TO 32u0 

00620 

224* 

3150 

IF(XT2.GT.VTIM(JrN> ) GO TO 3200 

00622 

225* 


K2« I)=0.0 

00623 

2*6* 

3200 

CONTINUE 

0 0 o 2 4 

227* 


TSTR(I*U=(K1(I)*IEXP(-Q(I)*XT1)-EXP(-Q(I)*XT2) ) ♦K2(I)*TIME(I) )* 

0 0624 

228* 


l ( lAbS(STRESSU) )**(EXPN(I)-1.) )*STRESS(I) >+TSTR(I-l»L> 

00625 

229* 

3400* 

CONTINUE 

00o26 

2j0* 

75 

T2<I*L)=TS+T1ME<I) 

® 0 0627 

231* 


TS=T2(IrL) 

>L00630 

23 2* 


GO TO 130 

O 00631 

233* 

300 

CREEP ( I » L ) = 0 • 0 

00t»J2 

234* 

130 

CONTINUE 

00634 

235 * 

131 

continue 

00636 

236* 


NR1TE < 6# 116) TITLE 

0 064 4 

237* 

116 

FORMAT(1H1#10A6> 

00645 

( 238* 


VkR1TE(6*141) 

00u47 

239* 

141 

FORMAT (1H0»4X»6HS1 HESS »12X»4HTEMP»13X.5HH0URS»13X»7HRUPTURE*8X»12H 

00647 

240* 

lTIME/RUPTUKE»7Xr8hSUM T/TR * 12X r 6HMARGIN/5X . 5H (PSl > » 13X. 3H<F ) * 13X » 7 

00647 

241* 

2HApPLltD»13X*4HTIi«it) 

00650 

2h2* 


FX=Q « 0 

00651 

243* 


CX=0.0 

00o32 

244 * 


DO 135 1=1 #NS 

Out 55 

245* 


FSl M(I )=FT(I)+FX 

00656 

♦DIAGNOSTIC* 

THE TEST FOK EQUALITY BETWEEN NON-INTEGERS MAY NOT BE MEANINGFUL. 

00o56 

246* 


IF (STRMAX.EQ.O.O) GO TO 275 

0iJo60 

247* 


DO 905 L=l»2 

0 0663 

24ft* 


DELSTW=STRMAX-TSTR(I»U 

0UO64 

249* 

905 

CONTINUE 

00666 

250* 


IF (DELSTR.LE.0.0) GO TO 280 

00670 

231* 

275 

CONTINUE 

0 06/1 

23 2* 


IF (NSIG.EO.O) GO TO 145 

006/3 

'* 253* 


XMAR ( I ) = ( 1 • O/FSUM 1 I ) ) -1 . Q 

00674 

234* 


IF (XMAR(I) .GE.O.U) 60 TO 145 

nuo 76 

235* 


toRITt (b. 150) XMAR(l) 

0U7(U 

2 36* 

150 

FORMAT (4tlH A NtGATIVE MARGIN HAS OCCURRED#MArG1N= »E15.8) 

00702 

237* 


N- 1 

O U703 

2 36 * 


00 10 lbb 

00/04 

239* 

145 

FXxFSl'Md ) 

0 0706 

2c0 « 


CX-Th.TR ( I .L) 
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00706 

2bl* 

00717 

2o2* 

00720 

2 63* 

00722 

2b4* 

00724 

2b5* 

00724 

2t»6* 

00724 

2b7* 

00/25 

260* 

00726 

2b 9* 

00732 

(270* 

00733 

271 + 

00742 

2 72 + 

00743 

2/3 + 

00/45 

274 + 

00745 

275 + 

00/45 

27.6* 

00746 

277* 

00752 

278 + 

0U76.1 

279* 

00763 

2b0* 

0U763 

2bl * 

00/63 

262* 

00764 

2o3* 

00770 

2b4* 

00777 

2o5* 

M 01o00 

2too* 

*01003 

2c7* 

•— >0 1003 

2b8* 

•“■O 1004 

2b 9 + 

01005 

29 0* 

01012 

291 + 

01012 

292+ 

01013 

293* 

01014 > 

294 + 

01015 

295* 

01020 

296* 

01021 

297* 

01U23 

290* 

01024 

299* 

01025 

3u0 + 

01030 

3U 1 * 

0103) 

3u2 + 

U10J3 

3u3* 

01037 

304* 

01037 

305* 

01040 

3ub* 

01041 

307 + 

01042 

308* 

01045 

309 + 

01047 

310* 

0 IU50 

311* 

01061 

312 + 

01053 

313* 

01053 

314 + 

. 'i i 1)64 

315 * 

0) Ot) 7 

316 + 

0 \ 

317 « 


WRITE (6*140) STRESS! I ) » TEMPI I). TIME (I) .TR i I j .FT { I ) .FSUMl I ) »XMAR( I ) 
140 FORMAT (lH0.7(E15«b*3X) ) 

135 CONTINUE 

WRITE (6*142) 

142 FORMAT ( 1H1 .86HTHE FOLLOWING VALUES’ ARE THE CREEP STRAIN VS TIME VA 
1LUES FOR THE TIME HARDENING THE0RY/5X.4HTIME.15X.12HCREEP STRAlN/5 
2X.4H(HR> » 16X »7H < IN/ IN ) ) 

TOSO.O 

WRITE (6* 143) T0.S1R0 

143 FCRMAT(lH0»E15.e.5A.El5.6> 

WRITE (6# 144) (T2(1.1).TSTR(I.1).I=1»NS) 

144 FORMAT (1H0.E15.6.SX.E15.6) 

WRITE (6. 146) 

146 FORMAT ( 1H1 .fibHYHt FOLLOWING VALUES ARE THE CREEP fRAlN VS TIME VA 
1LUES FOR THE strain HARDENING THE0RY/5X.4HTIME.15X.12HCREEP STRAIN 
2/5X >4H(HR ) » 16X . 7H ( IN/ IN ) ) 

WRITE (6. 143) TC.STRO 

WRITE <6. 144) <T2U»?)»TSTR(I»2).I=1.NS) 

WRITE (fa* 147 ) 

147 FORMAT (lHl.ftlHTHE HOLLOWING VALUES ARE THE CREEP STRAIN VS TIME VA 
1LUES For THE PAO-MARIN THE0RY/5X»4HTIME»15X»1?HCREEP STRAIN/5X»4H< 
2HR ) » 16X. 7H ( IN/IN ) ) 

WR1TEI6.143) TO .GIRO 

WRITE <6» 144) (T2(1.3)*TSTR(I»3) » 1=1 .NS) 

GO TO 165 

260 WRITE! 6.290 ) TSTRU.L) 

290 FORMAT ( 1H0.55HTOTAL STRAIN HAS EXCEEOED MAX STRAIN ALLOWABLE* STRA 
1IN=,E15.8) 

GO TO 165 

155 WRITE (6. 160) STRESS (N) »TEMP(N) »TIME(N> 

160 FORMAT (2HH RUPTURE OCCURRED AT STRESS=»E15,8.5X5HTEMP=»E15.8.5X»5H 
1TIME=.E15.8> 

165 CONTINUE 
TISUM=0.0 
DO 400 1 = 1 * NS 
TISUM=TISUM ♦ TIKL(I) 

400 CONTINUE 

XL1F£= TISUM/FSUM INS ) 

TMAXSTEMPtl) 

WRITE(6*425) XLIFE 

425 FORMAT ( 1H1 »5X» 14HLIFE EXPECTE0=*E15.8) 

IF (XLIFE. OT.L1REO) GO TO 455 
WR1 TE ( fo » 460 ) XLIFt.LIRtO 

460 FORMAT (1H0.17HTHE LIFE EXPECTED.E15.8t 31H IS LESS THAN THE LIFE RE 
1GUIRED.E15.8) 

GO TO 465 
455 CONTINUE 

00 410 1=1. NS 

IF ITFAX.GE.TEMPm ) GO TO 420 
TMAX=TEMH(1) 

420 CONTINUE 
410 CONTINUE 

COEFEM=(PCR(3)-PCKl?) )*1000./( ( TMAX+460 . ) * ( ALOGIO ( S1GCR (2 ) )-ALOG10 
1 (S1GCR (3) ) ) ) 

WRITE (<’»430) TMAX 

430 FORMAT ( 1h;Q.5X.9HMmX TEMU=E15.8) 

WRITE (n.4’35) CoEFNM 
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DATE 010373 

01063 318* 435 FORMAT ( 1H0*5X»2HM=»E15.8) 

01064 *DIAGN0STIC* THE Tt ST FOR EQUALITY BETWEEN NON-INTEGERS HAT NOT BE MEANINGFUL. 


01064 

319* 


IF(LIREQ.EU.O.O) GO TO 436 

01066 

320* 


FS= ( XL1FE/LIHEG > ** (1 . O/COEFFM ) 

01067 

321* 


WRITE (6.440) FS 

01072 

322* 

440 

FORMAT(1HO»5X»17H(-ACTOR OF SAFETYCE15.8) 

01073 

323* 

436 

CONTINUE 

01074 

324* 

465 

CONTINUE 

0 1075 

325* 


NP=NS+1 

01076 

,326* 


HO 800 1=1 »NP 

01101 

3^7* 


IF (l.GT.l) GO TO 805 

01103 

328* 


TX ( 1 . 1 ) =T0 

01104 

329* 


SY(l,l)=bTKO 

01105 

330* 


GO TO 60U 

0x106 

331* 

805 

continue 

01107 

332* 


TX ( I > 1 ) =T2 ( 1-1 . 1 ) 

oniu 

333* 


SY(I.1)=TSTR(I-1»1) 

01111 

3 34 * 

800 

CONTINUE 

01113 

335* 


CALL IDENT (9) 

01114 

336* 


CALL QUIK3V (1.35* t#COX » BCD Y » -NP .TX.SY) 

01115 

337* 


DAiA BCD Y /12*bH / 

01117 

338* 


DAI A PCDX / 6HTIMt <»6HHK> .10*6H / 

01121 

339* 


CALL APRNTV(0.-14»20»2UHCREEP STRAIN < IN/lN> * 0 . 652) 

U X 122 

3h0* 


CALL PRlNTv(21.21hT IME HARDENING THEORY.80. 1007) 

01123 

341* 


DO 910 1=1. NP 

ft) 01126 

3h2* 


IF(I.GT.l) GO TO 920 

i 0)130 

3h3* 


TX(1»?)=10 

£01131 

344* 


SY(l»2)=bTK0 

Oil 32 

34 5* 


GO TO 910 

0 1133 

346* 

920 

CONTINUE 

01134 

347* 


TX ( I .2 ) =12 ( 1-1 »2) 

01135 

348* 


SY( I.2)=TSTR(I-1.2) 

01136 

349* 

910 

CONTINUE 

0 114 0 ! 

330* 


CALL CUIK3V (-1 .3b.BC0X.BCDY .-NP.TX (1.2) .SY11.2)) 

01141 

351* 


DATA BCDY /12*6H / 

0.1 143 

352* 


DATA PCDX /6HTIVE 1.6HHR) »10*6H / 

01145 

353* 


CALL APRNTV(0*-14 » 20 . 20HCREEP STRAIN ( IN/IN) »0»652) 

01146 

3d4* 


CALL PRINTV(23»23HSTRA1N HARDENING THEORY » 60 . 1C07 ) 

01147 

355* 


DO 1010 1=1. NP 

01162 

366* 


IFU.GT.l) GO TO 1020 

01164 

357* 


TX ( 1 . 3 ) =TO 

01155 

358* 


SY(1.3)=STRO 

Oil 66 

359* 


GO TO 1010 

01167 

3oO* 

1020 

CONTINUE 

01160 

3d 1 * 


TX(I. 3) =12(1—1.3) 

01161 

3o2* 


SY(I.3)=TSTR(I-1 .0) 

01162 

3o3* 

1010 

CONTINUE 

0 1164 

3 o 4 * 


CALL GUIK3V(-1» 35 .bCDX .BCDY » -NP.TX (1.3) »SY(1.3) ) 

01165 

3o5- 


DATA BCDY / 1 2*6H / 

01167 

306'* 


DATA PCDX /6HT IMt ( »6HHR ) .10*6H / 

01171 

3o7 * 


CALL APRNTV ( 0.-14. 20. 20HCREEP STRAIN ( IN/IN) . 0.652) 

011/2 

3c6* 


CALL FRINTV(16»lbhHAC-MAhIN 1 HEOHY » PO . 1 007 ) 

0 i 173 

3u9^ 


CALL GUIn 3V (-1.24. oCDX .HCOY .-NP.TX ( 1 . 2 ) »bY(l»?) ) 

U 174 

3/0* 


DMA Kt'Y / 1 2*6H / 

H J t /Iy 

3 / 1 * 


D ft i A HLO* /(.hTIvl (fi.HMP) . 1 0*()H / 

0 i £.'00 

3/2* 


CALL APRNTV ( 0 » -1 4 . 20 . InHCPEF.F STRAIN ( IN/IN) . 0 , 652 ) 

i'l Jill 

3 / 5 * 


CALL f '.I6IV < U >2 r . i^LL'X m-.CLY . -NP.TX (l.l)»SY(t»l)> 


PAGE B 
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01202 

374* 

data 

BCDY /12*6H / 

0120** 

375* 

BA T A 

BCDX / 12*6H / 

01206 

376* 

CALL 

G)UIK3V(0'21'bCDX’bCDY'>NP'TX('l'3)Uf< 

01207 

377* 

DATA 

BCOY /12*6H / ■ T 

01211 

37B* 

data 

BCDX /12*6H / 

01213 

379* 

CALL 

F.NDJOB 

01214 

3b0* 

END 



END OF 

compilation: 

b DIAGNOSTICS. 


w 


DATE 010373 


PAGE 


i Cl»3)) 
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SFOR* lb OlK 

FOR 0£ : 10M-01'03/73-li:ii:34 (.0) 


SUt'MOuT XNl£ 6IR1 ENTRY POINT 000100 


STORAGE USED} CODE(l> 0001l6» DATA(O) .0000261 BLANK C0MM0N(2> 000000 


fcXTtRNAL REFERENCES (BLOCK. NAME) 
00u3 NERR3S 


STORAGE .ASSIGNMENT (BLOCK. TYPE. RELATIVE LOCATION. NAME) 

0001 000014 105G 0001 000024 20L 0001 000031 30L 

COUO K 000001 A2 0000 R 000005 UA 0000 R 000006 D1 

0000 000010 1NJP$ 0000 I 000003 11 0000 I 000002 12 


, OOlUl 

1* 

SUBROUTINE GIR1 (A.B.C.U.N) 

00103 

2 * 

DIMENSION C(50) .U(50) 

, 00104 

3* 

DO 10 1=1. N 

00107 

4* 

IF(A-C<I) >10.20.30 

00112 

5* 

10 CONTINUE 

00114 

6* 

20 B=D(I) 

0U115 

7* 

GO TO 40 

00116 

8* 

30 A2=C(I) 

00117 

9* 

12=1 

00120 i 

10* 

11=1-1 

00121 

11* 

A1=C (ID 

00122 

12* 

DA=(A1-A)/(A1-A2) 

00123 * 

13* 

D1=D( I 1 ) 

00124 

14* 

D2=D(I2) 

00126 

16* 

B=Dl*t»A* (D2-D1 ) 

00126 

16* 

40 RETURN 

00127 

17* 

ENU 


END OF COMPILATION 


NO DIAGNOSTICS 


DATE 010373 


PAGE 


10 


0001 000065 40L 0000 R 000004 A1 

0000 R 000007 D? 0000 I 000000 I 
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CFOR'lb IN|ER 

FOR OBlOA-Ul/03/73-li: 11736 t*0> 


DATE 010373 


RAGE 


11 


SUBROUTINE 1NTER1 ENTRY POINT 000102 

STORAGE USEO: C0DE(1> 000121l OATA(O) 000030) blank common ( 2 ) 000000 

( 

EXTERNAL REFERENCES (BLOCK * NAME) 

0003 NEKR3S 


STORAGE ASSIGNMENT (BLOCK* TYPE* RELATIVE LOCATION* NAME) 


COul 000014 10SG 0001 000024 20L 

0000 k oooool A2 0000 R 000005 L) A 

COUO 0J0010 INJPS 0000 I 000003 11 


00101 

1 * 


bd 00103 

2* 


i 001 Oil 

3+ 


£g0107 

4* 


00112 

5* 

10 

00114 

6* 

20 

00115 

7 * , 


00116 

e* 

30 

00117 

9* 


00120 1 

10* 


00121 

11* 


00122 

12* 


00123 

13* 


00124 

14* 


00125 

lb* 


00126 

16* 

40 

00127 

17* 



SUBROUTINE INTLR1(A*B*C*0*N> 
DIMENSION C(20) *0(20) 

00 10 1=1 *N 
IF (A-C (I)) 30 *20 *10 
CONTINUE 
B=D ( 1 ) 

GO TO 40 
A2=C ( I ) 

12=1 
11 = 1-1 
A1=C (ID 

DA=(a2-A)/(A2-A1) 

D1=0(I1) 

02 = 0 ( 12 ) 

B=l)2-0A* (02-01) 

RETURN 

END 


0001 000031 30L 

0000 R 000006 01 
0000 I 000002 12 


0001 000066 40L 

0000 R 000007 02 


tND OF COMPILATION: no DIAGNOSTICS. 


!»RRT *T 

FUR PUR 023A— 01/03-11 * 11 
CKttPAi,ALYSI*TPFS tLEMENT TABLE 


KiAMt 

VtRSION 

TYPE 

DATE 


time: 

SEO tt 

SI2E-PKE# TEXT 

(CYCLE 

WORD) 

mA I u 


FUR SYMb 

03 JmN 

73 

li:ii:25 

1 

QO 

5 

0 1 

r/A In 


RtLOCATAPLL 

03 JAN 

73 

li:ii:33 

2 

2 103 



G XU 


l-OR SYMb 

03 JAN 

73 

li:ii:34 

3 

3 

5 

0 1 

G I K 


HLLOCATAPLc. 

03 JAN 

73 

11 : 1 i:36 

4 

1 b 




0000 R 000004 A1 
0000 I 000000 I 


PSRKODE LOCATION 
1792 
1882 
1987 
1990 
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SECTION B.4 PROGRAM OUTPUT - SAMPLE PROBLEM 

OATE 010373 PAGE 26 

MANSON-SUCCOP PARAMETER FOR CREEP RUPTURE FOR RENE 41 BAR 

» lli0u0U00+06 • 16250000+02 .80000000+05 . 16900000+02 x/*io000o0o+05 .17750000+02 .40000000+05 .18700000+02 

4r 

.30000000+05 .19300000+02 .20000000+05 .20100000+02 .15000000+05 .20700000+02 .10000000+05 .21300000+02 

LARSuN-MLLER PARAMETER FOR .2 0/0 CRFEP»RENE 41 BAR 
AMOUNI OF, CREEP STRAINS .20000000-02 

. lii5u0U00+06 .36000000+02 .10000000+06 .36200000+02 .80000000+05 .37700000+02 .60000000+05 ,39550000+02 

• 4 l)0uUO00 + 05 .41850000+02 .20000000+05 .44750000+02 .10000000+05 .46800000 + 02 

MAX STRAIN ALLOWABLES, .15000000+00 
LIFE REuUiRED- .20000000+03 


a 

l 

i— * 

O' 
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DATE 010373 PAGE 27 

THE FOLLOWING VALUES ARE THE PRESCRIBED LOADING CYCLE 
STRESS(HSI) TEMP(F) TIME (HR) 

.250u0000+05 .16000000+04 .40000000+01 

.35u0u0u0+05 .15000000+04 .35000000+01 

,85u000o0+05 .13000000+04 .45000000+01 

• 10u000u0 + 06 .13000000 + 04 . 15000000+01*-^ 


to 
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OATE 010373 


PAGE 


28 


THE FuLLOfcING OATA ARE 

: the transient 

CREEP STRAIN VALUES FOR 

THE LOAD STRESS AND 

temperatupe values 


NV= 9 STRESS* 

*25000000405 

TEMP* .16000000*09 




.OOuOOOOO 

.00000000 

.10000000400 

-.10000000-09 

.10000000*01 

.5000u000-09 

.50o0U0U0+0l 

.25000000-03 

.60000000401 

.90000000-03 

.10000000*02 

.48000000-03 

.Uu0u000+02 

.52500000-03 

,12000000+'!2 

.55000000-03 

.13000000*02 

.57000000-03 

NVs 11 STRESS* 

.65000000405 

TEMPS .13000950409 




. OOoOuOoO 

.00000000 

.10000000+00 

.90000000-09 

.50000000+00 

.15000000-03 

.lOuOUOUO+Ol 

.25000000-03 

.30000000*01 

.52000000-03 

.50000000+01 

.75000000-03 

.6nuOoooo+oi 

.85000000-03 

.70000000*01 

.92000000-03 

.80000000*01 

.96000000-03 

.9OoOU0oO+01 

.99000000-03 

.10000000*02 

.10000000-02 



NV= lj STRESS* 

.35000000405 

TEMPS .15000000*09 




• UOuOuOuO 

.00000000 

.10000000*00 

.90000000-09 

.50000000+00 

.15000000-03 

. lOuOUOuO + Ol 

.25000000-03 

.20000000+01 

.99000000-03 

.90000000+01 

.80000000-03 

•55U0U0G0+01 

.10500000-02 

.75000000*01 

.13500000-02 

.90000000*01 

.15000000-02 

• 12u000u0+02 

.17000000-02 

.13000000*02 

.17500000-02 

.14000000+02 

.17700000-02 

• lb 00 u 0 u 0 + 02 

. 16000000-02 





NV= o STRESS* 

.10000000406 

TEMP= .13000000*09 




. OOuOuOOO 

.00000000 

.10000000+00 

.10000000-03 

•50000000+00 

.32500000-03 

. lOyOuGUO+Ol 

.56000000-03 

.15000000*01 

.79000000-03 

.17500000+01 

.6900u000-03 


W 


00 
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DATE 010373 


PAGE 


SAMPLE PROBLEM- TIME AND STRAIN HARGEN1NG.TRANSIENT ♦ 

STRtSS 
<PS1> 


.25000000+05 
•350O0UO0+O5 
.85000000+05 
• lUOUOoOO+06 


TEMP 

<F> 

.16000000+04 

.15000000+04 

*13000000+04 

•13000000+04 


hOURS 

APPLIED 

.40000000+01 

.35000000+01 

.45000000+01 

.15000000+01 


SO. 

RUPTURE* 

TIME? 

.23473314+03 

.60047882+03 

.48230762+03 

•16218102+03 


TIME/RUPTURE 

.16685219-01 

.58286818-02 

•93284041-02 

•92489241-02 


SUM T/TR 

.16685219-01 

.22513901-01 

.31842304-01 

.41091228-01 


MARGIN 

.58933286+02 

.43417003+02 

.30404762+02 

•23336094+02 



B-20 


I 


THE FuLLC*IN6 VALUES 
TiMt • 

<HRJ 

.OoOOOOOO 
•AOUoGUOO+Ol 
• 7&000.000+01 
« l20uGUO0+O2 
•13500000+02 


ARE the creep strain vs 

CREEP STRAIN 
(IN/IN) 

*00000000 

*20000000*03 

.75000001*03 

*10609137*02 

. 10797622*02 


TIME VALUES FOR THE 


TIME hardening THEORY 


OATE 010373 PAGE 30 
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DATE 010373 PAGE 31 


THE FOLLOWING VALUES ARE THE CREEP STRAIN VS TIKE VALUES FOR THE STRAIN HARDENING THEORY 
TIME CREEP STRAIN 

tHR) UN/IN) 


.ooooouoo 
»4u0u0009+01 
• 7;>0u0U00+01 
.Ii0o0000+02 
.13500000+02 


.00000000 

.20000000-03 

.04106667-03 

.10539404-02 

.18647889-02 


w 
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following values 
time 

<HR) 

are the creep strain vs time values for the pao-marin theory 

CREEP STRAIN 
(1N/IN) 

. OOOuOOOU 

.00000000 

.4uoooooo+oi 

.20000000-03 

• 7t)000o00+01 

.83764221-03 

• I20o0u00+02 

.10279986-02 

. liSu0U00+02 

.18388472-02 


OATE 010373 PAGE 32 


i 
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line HARDEN INli 1HE0RV 
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Fig.B-2 - Creep Strain, Time Hardening Theory 


B-23 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 




LMSC-HREC TR D306579 


STRAIN HARDENING THEORY 
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Fig.B-3 - Creep Strain, Strain Hardening Theory 
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mchiarin iheohy 



Fig.B-4 - Creep Strain, Pao-Marin Theory 
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